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listing  was  prepared  and  included  as  part  of  this  report.”  The  Task  2  effort  re¬ 
sulted  in  a  model  for  estimating  the  number  of  fragment/debris  missiles  per  unit 
of  ground  surface  area  which  exceed  specified  impact  energy  levels.  The  Task  3 
effort  resulted  in  a  scale  model  test  program  and  similitude  analysis  whicli  can 
be  used  to  generalize  the  test  results. 
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FOREWORD 


This  report  documents  the  results  of  an  analytical  program  performed  by  the 
Southwest  Research  Institute  (SwRl).  The  program  was  intended  to  provide  some 
of  the  necessary  methodologies  and  data  required  to  begin  implementation  of  th3 
Naval  Ordnance  Hazards  Analysis  and  Risk  Management  (NOHARM)  System  during  the 
1980-81  timeframe.  NOHARM  is  a  software  system  developed  by  the  Naval  Civil 
Engineering  Laboratory  (NCEL)  (W.  Keenan,  Program  Manager).  It  is  designed  to 
provide  information  required  to  assess  and  manage  the  risks  to  persbnne'  and 
property  exposed  to  hazards  associated  with  the  handling  of  naval  ordnance. 

This  effort  was  performed  by  Southwest  Research  Institute  (J.  C.  Hokanson, 
principal  investigator)  for  the  Naval  Sur  ace  Weapons  Center  (NSWC)  under 
Contract  Number  N60921-80-C0267.  Dr.  J.  M.  Ward  (NSWC)  was  the  contract 
monitor.  Dr.  L.  Huang  (NCEL)  monitored  the  NSWC  effort  performed  under  Contract 
N68305-80-WR00101. 


This  report  does  not  reflect  the  official  view  or  final  judgement  of  NSWC. 
Its  main  purpose  is  to  publish  the  review  of  technical  literature  and  the  SwRI 
approach  for  the  predictive  methodologies  and  master  test  plan  for  estiroatinjg 
debris  hazard  environments  from  reinforced  concrete  structures  subjected  to 
internal  explosion  loading. 

The  technical  work  was  completed  in  July  1981  and,  therefore,  does  not 
reflect  any  of  the  advances  in  fragment  and  debris  hazard  prediction 
methodologies  which  have  been  developed  since  July  1981. 
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CHAPTER  I 
INTRODUCTION 


This  report  presents  the  results  of  an  analytical  program  intended  to 
provide  necessary  methodologies  and  data  required  to  begin  implementation  of  the 
Naval  Ordnance  Hazards  Analysis  and  Risk  Management  (NOHARM)  system.  This 
effort  was  conducted  for  the  Naval  Surface  Weapons  Center  (NSWC),  White  Oak, 
under  Contract  Number  N60921-80-C-0267.  NOHARM  is  a  software  system  designed  to 
provide  information  required  to  assess  and  manage  the  risks  to  personnel  and 
property  exposed  to  hazards  associated  with  the  handling  of  naval  ordnance.  The 
primary  objective  of  the  NOHARM  system  is  to  identify  unsafe  (unacceptable  risk) 
conditions  in  ordnance  handling.  A  secondary  goal  is  to  determine  the  optimum 
strategy  for  mitigating  risks  once  they  are  discovered.  To  achieve  these 
objectives,  NOHARM  must  have  the  capability  to  estimate  the  explosive  hazards 
associated  with  ordnance  operations,  predict  the  human  and  economic  risks  from 
possible  explosions  and  fires,  assess  those  risks  to  identify  unsafe  conditions, 
and  assist  in  the  selection  of  optimum  techniques  for  risk  mitigation.  As  such, 
tlie  software  will  consist  of  three  distinct  modules:  the  Explosives  Hazards 
Model  (EHM),  the  Risk  Prediction  Model  (RPM),  and  the  Risk  Mitigation  Model 
(RMM).  This  report  is  concerned  only  with  the  EHM.  The  EHM  is  intended  to 
produce  several  primary  outputs  which  are  required  by  the  other  two  modules. 
These  outputs  are: 

o  The  probability  per  year  of  all  possible  yields  of  explosions  and 
fires  from  transactions  and  storage  of  naval  ordnance. 

o  The  maximum  credible  yield  from  such  accidents. 

o  The  identification  of  particular  ordnance  transactions  which  are 
major  causes  of  fire  and  explosion. 

o  The  estimation  of  the  blast,  fragment,  and  debris  environment 

associated  with  explosions  in  the  naval  ordnance  handling  system. 

This  report  is  concerned  only  with  the  estimation  of  the  fragment  and  debris 
environment  around  accidental  explosions  in  naval  facilities.  This  project 
consisted  of  three  tasks;  a  comprehensive  review  of  the  technical  literature, 
the  development  of  a  new  methodology  for  predicting  the  debris  environment  from 
internal  explosions  in  reinforced  concrete  structures,  and  the  development  of  a 
master  plan  to  derive  the  input  data  required  for  full  implementation  of  the  EHM 
of  NOHARM. 

The  objective  of  Task  1,  literature  review,  was  to  prepare  a  bibliography 
on  the  primary  fragment  hazard  from  accidental  detonations  of  Navy  weapons,  and 
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the  debris  hazard  from  accidental  exploaions  both  in  reinforced  concrete 
buildings  and  in  ships. 

The  review  of  the  literature  on  fragment/debris  characteristics  was 
conducted  with  the  aid  of  three  computerized  information  retrieval  systems.  To 
supplement  the  computerized  scan  of  the  literature,  the  report  files  at 
Southwest  Research  Institute  (SwRI)  were  examined.  J.  M.  Ward  at  NSWC  also 
contributed  many  references  to  pertinent  reports.  Finally,  we  examined  the 
proceedings  of  the  Department  of  Defense  Explosives  Safety  Board  (DDESB) 
meetings  for  the  past  12  years.  Through  this  exhaustive  scan  of  the  literature, 
102  useful  reports  and  papers  f’ere  identified.  A  bibliographic  listing  of  these 
reports  was  prepared  and  is  included  as  Appendix  A  of  this  report. 

The  objective  of  Task  2  was  to  develop  a  methodology  for  predicting  the 
fragment  and  debris  environments  from  accidental  explosions  within  reinforced 
concrete  structures.  The  model  developed  provides  a  technique  for  the 
estimation  of  the  number  of  fragments  (emanating  from  the  weapon)  and  the  number 
of  debris  missiles  (emanating  from  the  concrete  structure)  per  square  foot  of 
ground  surface  area  which  exceed  specified  energy  levels  on  impact.  The 
methodology  is  presented  in  a  form  suitable  for  eventual  conversion  to  a 
CO jputer  code.  However,  certain  critical  input  parameters  are  not  currently 
available,  so  development  of  the  computer  code  is  restricted  to  only  those 
situations  for  which  experimental  data  are  available. 

The  methodology  presented  parallels  a  similar  methodology  developed  by 
J.  M.  Ward  for  predicting  the  fragment  hazard  associated  with  the  accidental 
detonation  of  weapons,  specifically  a  pallet  load  of  Mk  82  bombs.  Kis 
methodology,  which  is  an  application  of  concepts  developed  by  F.  B.  Porzel  for 
the  Naval  Explosive  Safety  Improvement  Program  (NESIP),  was  developed  for  the 
special  case  of  an  explosion  in  the  open  field  providing  for  both  high  and  low 
trajectories  of  the  debris  or  fragments.  This  feature  is  important  in  the 
establishment  of  acceptable  hazard  arcs.  In  this  report  we  have  included  this 
concept  in  the  estimation  of  hazards  for  both  fragments  and  debris  missiles. 

Task  3  was  intended  to  develop  a  master  plan  which  could  be  used  by  the 
Government  to  provide  necessary  methodologies  and  data  required  to  implement  the 
NOHARM  system  for  fragment  and  debris  effects.  The  objective  of  this  task  was 
to  identify  gaps  in  the  current  state-of-the-art  of  predicting  the  fragment  and 
debris  environment  from  accidental  explosions  of  naval  ordnance. 

During  this  project  we  have  established  that  major  data  gaps  exist  in  the 
prediction  of  structural  failure  patterns,  in  the  estimation  of  internal  loading 
on  the  buildings,  and  in  the  prediction  of  initial  cr\ditions  used  to  calculate 
debris  trajectories.  Much  weapon  effectiveness  data  are  available  and  are 
useful  for  establishing  fragment  effects  (from  the  weapons  case).  However,,  not 
much  data  for  the  debris  effects  (from  the  surrounding  structures)  are  available 

Currently,  to  predict  the  debris  hazards  for  a  particular  situation,  tests 
must  be  conducted  to  provide  the  requisite  data.  This  report  presents  a  scale 
model  test  program  and  similitude  analysis  which  may  be  used  to  generalize  the 
test  results. 
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CHAPTER  2 

REVIEW  OF  EXISTING  DEBRIS  AND  FRAGMENT  PREDICTION  METHODOLOGIES 


Thi  technical  literature  was  reviewed  to  identify  existing  methodologies  to 
define  the  hazards  associated  with  the  accidental  detonation  of  weapons  within 
and  outside  of  different  classes  of  structures.  Of  particular  interest  was  the 
explosion  of  naval  weapons  within  reinforced  concrete  structures  and  ships.  In 
the  following  paragraphs,  the  methodologies  identified  during  this  effort  are 
described  for  bomb  fragments  and  for  debris  from  ships  and  reinforced  concrete. 
In  this  report  we  will  use  the  term  fragment  to  refer  to  pieces  of  the  weapon 
and  the  term  debris  to  refer  to  pieces  of  the  structure  emanating  from  the 
explosion. 

The  review  of  the  literature  on  fragment/debris  characteristics  was 
conducted  in  several  stages.  The  COMPENDEX  (Computerized  Engineering  Index)  and 
the  National  Technical  Information  Service  (NTIS)  data  bases  were  scanned  using 
the  Lockheed  DIALOG  information  retrieval  system.  Addditionally ,  a  search  was 
initiated  on  the  Defense  Technical  Information  Center  (DTIC)  data  base.  To 
supplement  the  computerized  scar  of  the  literature,  the  report  files  of  several 
people  at  SwRI  were  examined.  J.  M.  Ward  also  contributed  many  references  to 
pertinent  reports.  Finally,  we  examined  the  proceedings  of  the  DDESB  meetings 
for  the  past  12  years.  Through  this  exhaustive  scan  of  the  literature,  useful 
reports  and  papers  were  identified.  A  bibliographic  listing  of  these  reports  is 
contained  in  Appendix  A. 

BOMB  FRAGMENT  t ITHODOLOGIES 

Existing  methodologies  for  predicting  the  behavior  of  bomb  fragments  have 
been  reviewed.  This  research  resulted  in  identifying  methodologies  for  the 
determination  of  each  of  the  following  parameters:  mass,  weight,  number,  size, 
velocity,  range  and  trajectory,  distribution,  and  impact  probability.  A  summary 
of  the  methodologies  found  for  the  determination  of  each  fragment  parameter  is 
presented  in  the  follc'.'ing  paragraphs.  Tables  are  included  to  give  specific 
methodologies  and  the  source  from  which  they  were  taken. 

Fragment  Mass 

The  most  extensive  treatment  of  bomb  fragment  masses  is  a  statistical  study 
conducted  by  Hekker  and  Pasman.^  An  equation  derived  by  Weibull  is  presented 
which  predicts  the  probability  that  a  mass,  x,  is  greater  than  a  given  mass,  m, 
This  equation  is  used  by  the  authors  to  derive  the  mean  fragment  mass,  m.  This 
term,  m,  is  i  function  of  m„,  which  is  also  the  average  mass  but  is  calculated 
differently  by  researchers  as  shown  in  Table  1.  It  is  important  to  distinguish 
between  m  and  mg  and  also  to  realize  that  other  researchers  such  as  Mott  use 
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TABLE  1.  METHODOLOGIES  FOR  PREDICTION  OF  FRAGMENT  MAtS 


PARAMETER 

METHODOLOGY 

DESIGNATION 

REFERENCE 

Probability 

P<x>m)  •  exp(-(n/m^)^] 

Uelbull 

(1) 

a^,  X  arc  functions  of  aacerlal  and  geometry 


Mean  Maas  x  •  ■  (1) 

a-  J x<iP{X>x)  -  a^r(I+l/.\) 

X  -  0 

where:  r  •  Euler  Caoaa  Function 


Average  Haas 


where:  •  Initial  fragment  velocity 

t  •  wall  thickness 

d . ,  d  •  diameters 

1  0 

U  ■■  metal  to  high  explosive  ratio 

C  “  coefficient,  function  of 
metal  and  high  explosive 

.  3^5/6 

where:  B  •  constant 

■  Inside  diameter 
t  •  thickness 


1/2 

u 


where:  C/M  ■  mass  ratio 
A  ■  constant 


where:  t  •  thickness 


d  •  Inside  diameter 
C  •  material  constant 

9 


Gurney  and 
Saraousakls 

(1) 

Lindeaen 

(1) 

Mott 

(!' 

Weiss 

(1 

Mott  (6 


Gurney  and  (6) 

Saraousakls 


Magls  (2) 


4 


u  in  place  cf  mo  whereas  Hekker  and  Far  an  use  p  to  symbolize  Che  ratio  of 
metal-to-high  explosive  mass. 


Kach  of  the  possible  average  mass  a  luations  presented  by  Hekker  and  Pasman 
was  tested  statistically  using  the  F-teoC  with  experimental  data  derived  from 
tests  on  ring-type  cylindrical  shells.  Their  results  indicated  that  the 
formulae  of  Gurney-Sarmousakis  dna  Mott  aru  Che  most  applicable  under  the 
conditions  tested.  A  note  is  made  by  the  authors,  however,  that  completely 
different  results  could  be  obtained  with  an  increase  in  the  number  of 
experiments,  or  from  a  charge  in  the  X -value. 


The  literature  reviewed  for  this  project  Indicates  that  the  formulae  of 
Mott  and  Gurney-Sarmousakis  ..re  tie  present  day  standard  for  calculating  the 
average  fragment  mass.  There  does  exist  an  alternative  to  these  two  equations, 
the  Magis  equation,  which  was  used  by  Randers-Pehrson,^'  and  is  presented  in 
Table  1. 


Fragment  Weight 

Some  of  Che  articles  considered  concern  themselves  with  fragment  weight, 
rather  Chan  fragment  mass.  Obviously  this  represents  only  a  minor  change  in 
variable  definition  (by  a  factor  of  the  acceleration  of  gravity)..  The  U.S.  Army 
Tuc  inical  Manu; 1  5-1300^  and  other  Army  work  use  -an  equation  for  primary 
fr  'meat  weight  based  on  Mott's  equation.  The  Swedish  researchers,  Ericksson 
at  Arvidsson,^  also  use  a  va.'iation  of  tha  Mott  equation  to  determine  the 
c -al  weight  of  all  fragments  with  a  weight  m,  greater  than  a  specified  mj^. 

This  equation  differs  dramatic.ally  from  others  encountered  in. Chat  it  includes 
th-'  influence  of  the  percentage  of  carbon  present  in  the  steel  shell  casing  on 
thfi  fragment  weight.  Fragment  weight  methodologies  are  given  in  Table  Z. 

Nrrber  of  Fragments 

Various  approaches  available  for  calculating  the  number  of  fragments  with  a 
mess  greater  than  a  given  macs, are  presented  in  Table  3.  A  number  of  these  ate 
p  "esented  by  Hekker  and  Pasn.an^  and  are  incorporated  into  their  statistical 
analysis,  which  has  been  previously  discussed.  The  Mott  equation,  or  variations 
thereof,  is  by  far  Che  most  commonly  used  method  in  this  category  and  appears 
rany  times.  Johnson  and  Moseley^  have  presented  2D  and  3D  variations  on 
M.jtt's  formula ,, while  Sternberg"  defined  a  fragment  range  for  the  Mott 
equation  and  developed  other  equations  for  the  number  of  fragments,  which  ace 
V’.lid  where  Che  Mott  formula  is  not.  Krauklis  and  Bedford^  present  the  Payman 
method  of  analysis  along  with  modifications , on  Mott  and  Payman.  Porzel,^ 
after  discussing  Mott's  general  equation,  proposes  a  prediction  equation  for  the 
number  of  fragments,  based  on  the  idea  that  1 ragmentation  is  controlled  by 
fragment  length,  rather  than  by  mass  or  area. 

F-agment  Size 

The  s  andard  method  for  predicting  fragm'Tit  ^Ire  (area)  is  a  formula  vhich 
i.s  a  function  of  the  fragment's  mass  and  the  shape  factor,  or  ballistic 
density.^  This  equation,  given  below,  is  the  only  method  used  to  determine 
fragment  size. 
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TABLE  3.  METHODOLOGIES  FOR  PREDICTION  OF  NUMBER  OF  FRAGMENTS 
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where 


k  “  ballistic  density 
A  “  presented  area 
m  mass 
Fragment  Velocity 

An  extensive  number  of  reports  containing  fragment  velocity  prediction 
equations  were  reviewed  as  part  of  this  program  and  a  list  of  prediction 
equations  is  given  in  Table  4.  Tl.e  Gurney  equation,  which  is  applicable  to 
cylinders,  is  the  most  widely  accepted  method  for  predicting  fragment  initial 
velocity.  Several  variations  of  the  Gurney  equation  have  been  developed  and  are 
also  included  in  Table  4.  Randers-Pehrson^  modified  the  Gurney  equation  tn 
predict  the  nonsteady  state  velocity  of  a  fragment.  He  also  modified  the  Gurney 
equativ>n  to  account  for  an  open-ended  cylinder.  Variations  of  the  Gurney 
equation  also  exist  for  spheres  and  steel-core  cylinders,  and  are  shown  in  Table 
4.^  The  maximum  initial  velocity  for  a  fragment  from  a  cylinder,  sphere,  and 
a  steel-cored  cylinder  is  given  in  U.S.  Army  TM  5-1300.^  This  manual  also 
gives  an  equation  to  calculate  the  striking  velocity  of  a  fragment  on  a  target 
for  small  fragments  up  to  several  ounces  in  weight.  The  initial  fragment 
velocity  (based  on  experimental  photographic  velocity),  fragment  mass,  and  the 
drag  coefficient  can  also  be  determine  in  an  equation  found  in  anotherU.S. 

Army  publication.^*^  Lacher^^  presentt.  in  equation  for  fragment  velocity  as 
a  function  of  drag  factor  and  range. 

The  only  equation  found  to  predict  secondary  fragment  velocities  appeared 
in  a  U.S.  Army  Corps  of  Engineers  publication.^^  Shape  factors  are  given  for 
different  geometries,  and  ranges  of  applicability  are  also  shown. 

Fragment  Range  and  Trajectory 

The  methodologies  encountered  for  the  prediction  of  fragment  range  and 
trajectory  are  given  in  Table  5.  Zaker^^  solves  the  equation  of  motion  for  a 
fragment  by  separating  the  problem  into  a  gravity-free  solution  and  a  perturbed 
solution.  Schreyer  and  Romesberg^^  give  equations  for  calculating  the 
accelerations  in  the  x  and  y' direction  for  two  sets  of  conditions:  first, 
assuming  no  fragment  lift  and  second,  assuming  both  fragment  lift  and  drag.  The 
latter  is  incorporated  into  the  FRISB  computer  code,  which  is  used  at  SwRI  and 
ref.,rf3u  later  in  this  report. 

The  ir.tial  trajectory  angle  can  be  found  using  the  Taylor  equation  for 
steady  state  conditions  which  is  referenced  in  numerous  reports.  Randers- 
Pehrson^  has  modified  the  steady  state  Taylor  equation  to  account  for 
nonsteady  state  conditions. 

Various  authors  give  relationships  between  fragment  range  and  velocity 
which  are  shown  in  Table  5,^^  along  with  validity  criterion  and  any 
assumptions  made.  Fugelson  and  Rathmann^^  present  a  range  equation  which  is  a 
function  of  the  initial  launch  angle  and  the  ratio  of  the  terminal  free-fall 
velocity  of  the  fragment  to  its  initial  velocity.  A  graphical  method  for 
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predicting  fragment  range  is  contained  in  NASA  Contractor  Report  3023,^^  and 
is  shovm  in  Figure  1.  This  set  of  curves  was  developed  by  performing  a  model 
analysis  to  generate  dimensionless  parameters  which  describe  ranges  for  selected 
cases,  and  then  the  results  were  plotted  to  form  the  curves.  It  should  be  noted 
that,  in  generating  these  curves,  sev  ral  initial  trajectory  angles  were  used  in 
the  analysis  to  obtain  the  maximum  range  for  the  respective  fragments. 

Fragment  Distribution 

Fragment  distribution  methodologies  reviewed  are  given  in  Table  6. 

Klein, assuming  a  Mott  analysis,  derives  a  formula  for  the  total  number  of 
fragments  per  unit  solid  angle.  An  expression  for  angular  fragment  distribution 
is  given  by  Johnson  and  Moseley^  while  Fugelso,  et  al.,^^  derive  a  formula 
for  the  fragment  density  in  terms  of  the  range  and  azimuth.  An  expression  for 
zonal  fragment  density  is  used  by  the  U.S.  Army  in  their  publication.  For 
large  range  fragment  densities,  Fugelso,  et  al.,^  present  a  formula  which  is 
a  function  of  the  angle  of  incidence  at  impact  and  the  ground  range. 

Fragment  Impact  Probability 

The  literature  search  conducted  for  this  program  resulted  in  only  one 
formula  for  predicting  the  probability  that  a  fragment  would  impact  a  particular 
target.^®  This  formula,  given  below,  appeared  in  numerous  publications  in 
different  forms. 


P  =  1  -  exp(-qAf) 


where 


A'p  *  target  area 
q  =  areal  fragment  density 


DEBRIS  FROM  EXPLOSIONS  IN  NAVY  SHIPS 


The  literature  search  indicates  that  methodologies  for  predicting  debris 
characteristics  from  explosions  inside  Navy  ships  and  submarines  are  very 
limited.  No  methodologies  were  found  to  predict  the  number  of  debris  fragments 
which  would  result  from  an  explosion  in  a  ship  or  submarine.  Some  model  tests 
have  been  performed  at  NSWC  where  debris  data  bases  were  created  which  include 
debris  of  particular  masses  and  sizes  recovered. To  date,  the  data  bases 
have  not  been  statistically  analyzed  to  establish  prediction  equations  or 
schemes.  The  only  report  located  which  presents  analytical  predictions  of 
debris  characteristics  for  Navy  ships  or  submarines  is  an  SwRI  report  on  debris 
hazards  from  an  explosion  in  a  torpedo  tender  workshop. However,  the 
breakup  pattern  assumed  resembles  a  pressure’  vessel  explosion  and  this 
assumption  must  be  considered  before  applying  the  methods  described  in  the 
report.  Tests  done  since  publication  of  that  work  show  that  explosions  in 
torpedo  workshops  usually  involve  detonation  of  several  warheads  which  result  in 
-  different  breakup  than  assumed  in  the  SwRI  report.  It  is  still  feasible  that 
gaseous  explosion  products  could  build  up  sufficient  pressure  in  a  ship  to  cause 
an  explosion  similar  to  a  pressure  vessel  rupture,  but  one  needs  to  analyze  what 
type  of  explosion  is  possible  before  predicting  debris  characteristics  using  the 
methods  contained  in  the  SwRI  report  (not  verified  by  test).  Although  th«» 
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FIGURE  1.  SCALED  CURVES  FOR  FRAGMENT /DEBRIS.  RANGE  PREDICTION 
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TABLE  6.  METHODOLOGIES  FOR  PREDICTION  OF  FRAGMENT  DISTRIBUTION 
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applications  are  limited,  the  methodologies  will  be  presented  since  they  are  the 
only  ones  revealed  in  the  literature  search. 

Debris  Velocity  , 

The  velocity  of  debris  from  a  bulkhead  or  deck  in  a  tender  workshop  is 
determined  in  a  stepwise  manner  consisting  of  three  basic  phases.  The  first 
phase  consists  of  the  initial  impulsive  loading  from  the  blast  wave.  The  second 
nhase  involves  the  expansion  of  the  gases  produced  by  the  products  of  the 
explosion.  The  final  phase  occurs  as  the  debris  is  further  accelerated  by  the 
venting  of  the  gases  around  it  as  it  leaves  the  ship  and  begins  its  trajectory. 

The  extent  of  damage  and  a  breakup  pattern  are  established  for  a  680  kg 
(1500  lb)  charge  weight  by  examining  the  loading  on  individual  members  in  the 
bulkheads  and  decks  immediately  surrounding  the  charge.  The  report^^ 
concludes  that  these  bulkheads  and  decks  will  either  remain  intact  or  break  into 
large  piec<?.s  with  the  explosion,  of  680  kg  (1500  lb)  of  high  explosive  in  the 
tender  workshop.  Therefore,  the  additional  momentum  obtained  by  the  structural 
elements  due  to  gas  expansion  of  the  combustion  products  is  maximized  since  the 
gas  is  assumed  not  to  escape  between  debris  fragments.  The  average  reflected 
impulse  imparted  to  an  entire  exposed  bulkhead/deck  is  determined  at  the 
location  of  the  average  distance  of  the  bulkhead/deck  from  the  charge.  The 
procedure  to  determine  the  initial  velocity  of  debris  once  the  impulse  has  been 
determined  is  presented  in  Table  7.  The  first  step  is  the  simple  calculation  of 
velocity  due  to  the  impulsive  loading  on  the  structure.  Kinetic  energy  is  then 
calculated  in  the  second  step.  The  strain  energies  (S.E.)  are  determined  for 
all  structural  elements  in  the  exposed  bulkhead/deck  individually,  added 
together,  and  subtracted  from  the  kinetic  energy.  The  adjusted  Initial 
velocity,  Vq,  is  then  determined  as  shown  in  thi  second  step  in  Table  7.  The 
third  step  in  the  procedure  is  the  determination  of  the  increase  in  velocity  due 
to  the  expansion  of  the  gaseous  explosion  products.  A  small  computer  program, 
GASEX,  is  used  to  solve  simultaneously  the  equations  of  motion  and  the  pressure 
equations  for  all  bulkheads/decks  surrounding  the  tender  workshop  where  the 
charge  is  located.  The  sequence  followed  by  the  computer  code  is  outlined  in 
Table  8. 

As  the  final  step  in  the  velocity  prediction  method,  the  computer  code 
EvXCYL  is  used  to  obtain  the  increase  in  velocity  due  to  venting  of  the  gaseous 
explosion  products  around  the  debris  as  they  begin  their  trajectories.  This 
code  is  a  modification  of  a  program  developed  by  SwRI  for  NASA^^  to  determine 
the  velocity  of  fragments  from  a  bursting  cylindrical  pressure  vessel.  A 
complete  description  of  the  program  can  be  found  in  Reference  22.  The  initial 
velocity  determined  form  the  blast  loading,  gas  expansion,  and  gas  venting  is 
next  used  to  calculate  the  range  that  the  debris  travels  from  the  ship. 

Debris  Range 

Maximum  debris  ranges  for  a  specific  initial  velocity  and  mass  (with  a 
prescribed  area)  are  determined  by  assuming  appropriate  launch  angles  as  input 
to  the  computer  code  FRISB  developed  by  Baker,  et  al.^"^  Several  probab.e 
breakup  patterns  are  considered  in  the  analytical  report. The  arena  data 
described  in  References  20  and  21  can  be  used  to  compare  with  analytical 
predictions  so  they  can  be  supported  or  modified. 
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TABLE  7.  DETERMINATION  OF  DEBRIS  INITIAL  VELOCITY 

Stepwise  Method:, 

1.  Impulsive  loading 

V  ■  — — 

M 

where  1  “  reflected  specific  impulse 
A  “  area  of  bulkhead/deck 

M  *  mass  of  bulkhead/deck. 

2.  Since  kinetic  energy  of  each  structural 
element  is 

K.E.  =  1/2  MV^, 

use  adjusted  initial  velocity 


where  S.E.  **  strain  energy. 

3*  Computer  code  GASEX  to  obtain  increase  in 
velocity  due  to  expansion  of  the  gaseous 
explosion  products. 

4.  Computer  code  BOXCYL  to  obtain  increase  in 
velocity  due  to  venting  of  gaseous  explosion 
products. 
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TABLE  8.  "GASEX"  COMPUTER  CODE  METHODOLOGY 

1.  Input  the  area,  mass,  initial  velocity,  and  initial 
position  from  the  center  of  the  explosion  of  each 
of  the  six  bulkheads/decks. 

2.  Input  the  initial  volume  and  pressure  p  of  the  torpedo 
workshop. 

3.  Input  atmospheric  pressure  p^  which  is  constant,  a 

start  time  of  zero,  arid  the  time  Increment  At  between 

* 

calculations. 

4.  Calculate  the  position  of  each  bulkhead/deck  at  txme 
(t  +  At)  from 

A  (At)^ 

(t  +  At)  =  [p(t)  -  p^)  - 

Vj^(t)  At  +  Xj^(c) 

[Note  that  At  is  assumed  to  be  so  small  that  p(t)  is 
essentially  constant  during  the, time  interval  between 
(t)  and  (t  +  At) ] . 

5.  Calculate  the  velocity  of  each  bulkhead/deck  at  time 
(t  +  At)  from 

A.  (At) 

v^  (t  +  At)  =  [p(t)  -  p^]  — — -  +  v^(t) 

i 

[Same  note  as  in  4  above]. 

6.  Calculate  a  new  volume 'V(t  +  At)  from 

6 

V  (t  +  At)  =  Z  [x.(t  +  At)  -  x.(t)]M.  +  V(t) 

1=1  ^  ^  ‘ 
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TABLE  8.  (CONT.) 

7.  Calculate  the  new  Internal  pressure  at  time 
(t  +  At)  from 

p  (t  +  At)  -  (40) 

V(t  +  At)"* 

8.  If  p  (t  +  At)  <  p^  then  stop  the  calculations. 

9.  If  the  longitudinal  bulkheads  have  reached  the  position 
of  the  shell,  then  stop  the  calculations. 

10.  Set  t  ■  t  +  At,  Xj^(t)  ■  (t  +  At),  and  Vj,(t)  ■ 

(t  +  At). 

11.  Go  to  4  above  and  continue. 
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Debris  Maas  and  Size 


No  reliable  prediction  methodologies  have  been  located  for  determining 
debris  mass  or  size.  In  an  NSWC  report^^  an  experimental  average  mass  was 
determined  in  model  tests  of  explosions  in  nuclear  attack  submarines  to 
establish  safe  handling  arcs  around  the  submarines  during  pierside  topping  off 
operations.  All  that  is  presented,  however,  are  experimental  data  used  to 
determine  the  arc  distances. 

Debris  Distribution 

All  reports  studied  are  concerned  with  the  quantity-distance  criteria  of 
determining  hazardous  debris  patterns  by  establishing  the  numbers  of 
debris/fragments  of  kinetic  energy  equal  to  or  greater  than  80J  (58  ft-lb)  per 
55.7  (600  ft^).  The  two  experimental  papers^®* ^  present  plots  of 

areal  density,  with  Reference  20  plotting  the  distribution  by  density  per  55.7 
(600  ft^).  Kulesz,  et  al. presents  debris  density  using  the  same 
criteria  for  several  breakup  scenarios.  The  density  is  determined  following  the 
procedure  in  Table  9.  It  should  be  noted  here  that  this  procedure  was  based  on 
specific  assumptions  about  the  breakup  of  a  portion  of  a  particular  ship.  The 
angle  6  is  governed  by  the  extent  of  damage  along  a  ship's  length. 

DEBRIS  METHODOLOGIES  FOR  REINFORCED  CONCRETE  STRUCTURES 

Reinforced  concrete  (R/C)  structures  are  widely  used  as  high  explosive  (HE) 
storage  magazines,  HE  manufacturing  facilities,  test  facilities,  and  other 
structures.  Many  structures  are  designed  to  contain*  an  accidental  detonation 
of  its  contents  using  design  methods  such  as  those  specified  in  TM5-1300. 

For  large  quantities  of  HE  it  may  be  impractical  to  design  for  containment, 
particularly  for  operations  which  are  considered  low  risk  such  as  magazine 
storage,  aircraft  shelters,  packaging  buildings,  and  other  areas  where  no 
machining,  pressing,  or  other  potentially  hazardous  activities  with  HE  are  being 
performed.  If  there  exist  large  distances  between  a  building  containing  HE 
(with  operations  at  any  risk  level)  and  other  occupied  areas  then  the  structure 
may  not  be  designed  to  contain  an  accident.  The  definition  of  what  constitutes 
a  dangerous  situation  for  debris  is  not  universally  determined;  however,  the 
DDESB  has  set  a  standard  of  not  more  than  one  fragment  per  55.7  m^  (600  ft  ) 
with  an  impact  energy  of  80J  (58  ft-lb)  or  more. 

The  problem  of  determining  expected  debris  hazard  to  other  occupied  areas 
♦'corn  an  explosion  in  an  R/C  building  has  been  approached  using  analysis  of  test 
ad  accident  data  to  obtain  statistical  fits  for  debris  parameters  (initial 
velocity,  range,  distribution,  etc.)  and  empirical  methods  which  apply  equations 
o'  motion  to  calculate  debris  parameters.  Methodologies  for  debris  from  R/C 
sttactures  include  the  following  categories: 


♦Containment — design  of  a  structure  to  withstand  an  explosion  and  mitigate 
debris  and  blast  hazards  to  other  occupied  areas.  This  may  include  permanent 
deformation  of  the  wall  without  collapse. 
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TABLE  9.  PROCEDURE  TO  DETERMINE  THE  HAZARDOUS  DEBRIS  DENSITY 
1.  Arc  Length  ■  R9 

where  R  Is  the  radius  of  a  sector  out  from  the  charge. 

R  Is  actually  the  sum  of  the  range  of  a  fragment  and  the  distance  from 
the  center  of  the  explosion  to  the  origin  of  the  fragment  (like  the 
shell  or  bulkhead  on  a  ship). 


55.7m2 


fragment  (^1)  ®  (width  of  debris)/  "sj  55.7. 


4.  Effective  number  of  55.7  ra  squares  covered  =• 

2 

N^  =  (number  of  debris) (number  of  55.7  m  squares 
covered  by  one  debris  fragment). 


5.  Debrls/55.7 


No.  of  55.7  m  squares  along  arc  length. 
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9  Debris  Size  or  Mass  and  Quantity  o  Distribution 

o  Initial  Velocity  o  Final  Velocity 

o  Launch  Angle  o  Drag  Effects 

o  Range 

This  section  of  the  report  will  review  mechodologies  for  each  of  the  previously 
listed  categories  individually  and  describes  one  report  which  combined 
methodologies  in  several  categories  for  a  complete  analysis  to  predict  the 
debris  hazard  from  an  explosion  inside  an  R/C  structure.  Also,  a  list  of 
available  debris  data  bases  will  be  given  and  discussed. 

Debris  Size,  Mass,  and  Quantity 

Data  bases  from  scaled  tests  and  accident  investigations  indicate  that 
debris  size  from  an  R/C  structure  can  range  from  very  small  pieces  up  to  entire 
parts  of  a  building  such  as  a  wall  or  roof.  This  wide  range  of  debris  sizes  as 
a  percent  of  the  total  building  area  lends  itself  to  a  statistical  fit.'  Kulesz, 
et  al.,^^  conducted  a  search  of  DDESB  accident  reports  containing  debris  data 
and  performed  a  statistical  analysis  of  debris  weight.  The  accident  data 
included  data  from  seven  accidents,  and  one  result  of  the  analysis  is  a  plot  of 
cumulative  probability  distribution  of  debris  weight  for  three  energy  levels 
(see  Figure  2).  All  buildings  in  the  data  base,  except  one,  were  primarily 
reinforced  concrete.  This  review  of  methodologies  could  find  no  correlation 
between  debris  weight  and  range  for  a  given  energy  level  and  it  is  suggested 
that  debris  weight  could  be  assumed  log  normally  distributed  within  a  given 
debris  range.  Figure  2  shows  that  approximately  30  percent  of  the  debris  is 
0.454  kg  (1.0  lb)  or  less  and  approximately  75  percent  of  the  debris  is  4.54  kg 
(10  lb)  or  less  for  the  energy  levels  indicated.  A  4.54  kg  (10  lb)  cube  of 
concrete  is  approximately  a  12.7  cm^  (5.0  in.)  cube,  hence  75  percent  of  the 
debris  will  be  relatively  small.  From  Figure  2,  debris  weighing  22.7  kg  (50  lb) 
and  up  comprises  only  approximately  10  percent  of  the  total  debris.  The  report 
noted  that  the  three  curves  on  Figure  2  aie  nearly  parallel,  'nd  standard 
deviations  are  almost  equal  for  all  log  normal  distributions.  This  could 
indicate  that  it  may  be  possible  to  derive  a  scale  factor  from  the  energy  ratios 
and  the  magnitudes  which  are  related  to  the  mean  of  the  distribution. 

Ahlers^^  studied  an  accidental  explosion  at  the  Pantex  Ordnance  Plant  in 
1960  and  described  debris  size  distribution.  Although  this  report  does  not  give 
a  debris  hazard  methodology,  Ahlers  made  several  observations  from  the  Pantex 
data  base  including: 

o  Larger  debris  did  not  have  as  great  a  range  as  smaller  debris. 

o  At  all  distances,  less  than  6  percent  of  the  total  debris  was 

greater  than  1.36  kg  (3.0  lb).  At  ranges  larger  than  366  m 
(1200  ft),  no  debris  above  0.45  kg  (1.0  lb)  was  found;  at 
intermediate  ranges  of  183  to  366m  (600  to  1200  ft),  the  percentage 
of  debris  above  1.36  kg  (3.0  lb)  is  greater  than  that  at  distances 
shorter  than  183  ra  (600  ft). 
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o  The  largest  concrete  fragment  weighed  1814  kg  (4000  lb)  and  was 
found  approximately  43.7  m  (150  ft)  from  explosion  center,  all 
other  debris  was  approximately  90.7  kg  (200  lb)  or  less. 

o  There  were  no  concrete  fragments  above  1.3  kg  (3.0  lb)  found  at 
ranges  greater  than  274  (900  ft),  which  could  suggest  an  optimum 
design  point  for  debris  hazards,  provided  this  is  characteristic  of 
explosions  in  R/C  structures  in  general. 

Vargas,  Hokanson,  and  Rindner^^  conducted  a  study  to  determine 
fragmentation  characteristics  of  reinforced  concrete  and  masonry  dividing  walls 
subjected  to  close-in  blast  effects.  The  study  included  a  model  analysis  and  a 
series  of  model  tests  on  1/6  scale  cantilevered  and  three-side  supported  R/C 
wallti  and  full-scale  concrete  block  masonry  walls.  An  extensive  amount  of  data 
were  collected.  A  model  analysis  was  conducted,  and  th^  results  guided  the  data 
presentation.  The  conclusions  drawn  from  this  useful  experimental  effort 
included: 

o  Debris  emanating  from  the  interior  of  the  panel  comprises  40 

percent  of  the  number  of  concrete  fragments  produced  in  any  test. 
Debris  originating  from  the  outside  face  (facing  the  recovery  area) 
comprises  another  40  percent  of  the  concrete  fragments.  The 
remaining  20  percent  of  the  concrete  fragments  are  produced  from 
the  :^vCeptor  (charge  side)  of  the  panel. 

o  Debris  produced  in  these  kinds  of  experiments  can  be  classified  as 
either  "chunky”  or  "pancake"  in  shape.  The  average  range  of 
"chunky"  concrete  fragments  is  generally  20  to  50  percent  greater 
than  the  average  range  of  the  "pancake"  concrete  fragments  on  a 
given  test. 

Based  on  the  statistical  distribution  of  debris  range,  mass,  and  velocity, 
the  following  conclusions  can  be  drawn: 

o  Mass  and  range  distributions  in  the  format  of  Mott  Distributions 
for  arena  fragmentation  tests  were  prepared.  The  resulting 
distributions  for  debris  range  and  mass  are  qualitatively  similar, 
and  similar  observations  were  drawn.  If  all  other  parameters  are 
held  constant,  more  concrete  fragments  at  each  mass  level  and  more 
concrete  fragments  at  each  range  level  are  produced  when: 

a.  the  total  impulse  applied  to  the  panel  is  increased, 

b.  the  panel  compressive  strength  is  decreased, 

c.  the  reinforcement  spacing  is  increased,  or 

d.  the  number  of  supporting  edges  is  increased. 
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Initial  Velocity 

The  velocity  of  individual  concrete  fragments  due  to  shock  loading  varies 
and  depends  on:  (1)  the  magnitude  of  Che  excess  impulse  defined  as  Che  blast 
impulse  minus  the  flexural  impulse  capacity  of  the  element  (area  under  Che 
resistance-time  curve),  (2)  Che  mass  of  the  debris,  (3)  Che  location  of  Che 
concrete  fragment  prior  to  collapse,  (4)  Che  interaction  between  the  debris 
during  their  flight,  add  (5)  Che  strength  and  time  history  of  the  compressive 
stress  wave  transmitted  through  the  R/C  wall  or  roof  as  the  blast  wave  is 
reflected.  Although  the  velocities  of  individual  concrete  fragments  differ,  the 
average  translational  velocity,  V^,  of  the  debris  after  complete  failure  can 
be  approximated  from  Che  excess  impulse,  ig,  and  the  momentum  of  Che  wall  or 
roof  after  collapse.  The  equation  below^  provides  a  means  of  estimating  the 
debris  velocities  from  the  blast  impulse  and  a  knowledge  of  the  R/C  geometry. 


where 


(3) 


“  applied  unit  blast  impulse 

Ph  =  reinforcement  ratio  in  the  horizontal  direction 

dg  “  distance  between  the  centroids  of  the  compression  and  tension 
reinforcement 

fds  ®  dynamic  design  stress  for  the  reinforcement 
H  “  span  height 

Vi  =  maximum  velocity  of  the  post-failure  debris 


Cy  “  impulse  coefficient  ; 

Cf  »  post-failure  debris  coefficient 

A  simpler  expression  widely  used  is  a  momentum  balance  that  ignores  strain 
energy  and  equates  applied  specific  impulse  times  area  with. kinetic  energy  as: 


V. 

1 


i  A 
m 


(4) 


where: 

A  “  area  of  applied  impulse 


i  =  applied  specific  impulse 

Vargas,  et  al.,^^  discussed  earlier,  present  plots  of  the  largest  debris 
velocity  as  a  function  of  impulse  factor.  The  report  concludes: 
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o  There  existed  a  wide  range  of  velocities  and  launch  angles  in  every 
test.  The  majority  of  the  debris  particles  were  observed  leaving 
the  wall  in  a  direction  perpendicular  to  the  original  surface,  and 
moving  with  a  velocity  at  or  near  the  maximum  observed  velocity. 
Late  in  the  event,  other  debris  particles  were  observed  which 
traveled  at  off-normal  trajectories  and  at  somewhat  lower 
velocities. 


o  The  largest  velocity  appeared  to  be  independent  of  the  reinforcing 
bar  (rebar)  spacing  but  dependent  on  total  impulse,  effective  wall 
thickness,  and  the  restraint  conditions  (the  panels  supported  on 
three  sides  had  higher  velocities  than  those  supported  on  one  edge, 
all  else  equal).  Some  effect  of  concrete  compressive  strength 
below  27.6  MPa  (4000  psi)  was  observed. 

o  Debris  velocities  for  masonry  walls  were  lower  than,  those  for  R/C 
walls. 

Debris  generated  from  completely  or  partially  enclosed  R/C  structures  which 
collapse  due  to  shock  loading  can  undergo  acceleration  from  the  quasi-static  or 
gas  phase  of  the  blast  loading  as  well  as  the  shock  phase.  To  account  for  this 
acceleration  and  its  contribution  to  debris  velocity,  Kulesz,  et  al.,^^  used  a 
computer  program  called  GASEX  to  calculate  the  acceleration  due  to  gas  phase  for 
an  explosion  in  a  torpedo  tender.  Typically  each  wall  and  roof  are  treated  as  a 
panel  with  initial  velocity  due  to  shock  loading.  The  quasi-static  pressure  is 
entered  as  initial  pressure.  In  a  time  step  fashion  the  panels  are  allowed  to 
move,  a  new  volume  and  internal  pressure  are  calculated,  and  time  step  repeated 
until  internal  pressure  is  atmospheric.  The  GASEX  program  is  explained  in  more 
detail  in  Table  8.  Debris  originating  from  a  particular  wall  or  roof  will  have 
the  same  velocity  as  that  for  the  wall  or  roof  panel  used  in  GASEX.  Similar 
programs  called  BOXCYL  and  CYLIN  can  be  used  which  take  into  account  volume 
change  and  pressure  venting  relief  between  concrete  fragments. 23  Debris  can 
be  generated  from  completely  or  partially  enclosed  R/C  structures  which  do  not 
collapse  due  to  shock  loading  but  do  collapse  due  to  the  quasi-stetic  loading. 

To  account  for  this,  the  structure  is  treated  as  a  pressure  vessel  which 
ruptures  and  breaks  into  panels  that  are  accelerated  due  to  the  expanding  gas. 
Reference  27  used  this  procedure  in  an  R/C  aircraft  shelter  by  running  the 
program  called  CYLIN. 


Launch  Angle 

The  angle  between  the  initial  velocity  vector  of  a  concrete  fragment  and 
the  ground  is  called  the  launch  angle  of  the  fragment.  Although  debris  from  an 
R/C  surface  may  take  on  a  wide  range  of  launch  angles,  the  predominate  angles 
are  approximately  normal  to  the  surface.  This  conclusion  was  suggested  by 
Vargas,  et  al.,2^  where  the  predominate  launch  angle  of  the  higher  velocity 
debris  was  normal  to  the  test  panel  surface.  Merz^S  reports  that  debris 
dispersal  tends  to  be  normal  to  the  building  surface.  Moseley  and  Whitney  ^ 
performed  an  analysis  on  an  R/C  Norwegian  aircraft  shelter  while  concurrent 
model  tests  (1/20  and  1/100  scale)  were  being  performed.  The  analysis  used 
launch  angles  for  debris  in  the  analysis  which  were  close  to  perpendicular  to 
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each  surface.  These  launch  angles  were  obtained  through  review  of  films  of  the 
model  tests  which  indicated  a  wide  range  of  launch  angles  from  each  surface  but 
a  predominance  normal  to  the  panel. 

Debris  Range 

Kulesz,  et  al.,^^  used  similitude  theory  to  organize  a  data  base  and 
performed  a  statistical  analysis  on  scaled  debris  range.  The  report  includes  a 
plot  of  cumulative  probability  versus  range,  R,  and  also  nondimens ional  range, 

R,  (range  divided  by  square  root  of  average  presented  debris  area,  A)  for  three 
energy  levels  which  are  reproduced  here  as  Figures  3  and  4.  These  curves  allow 
one  to  predict  the  percentage  of  concrete  from  an  R/C  structure  that  will  travel 
a  certain  distance.  The  report  noted  that  the  three  curves  on  Figures  3  and  4 
are  nearly  parallel.  This  could  indicate  that  it  may  be  possible  to  derive  a 
scale  factor  from  the  energy  ratios  and  magnitude  which  is  related  to  the  mean 
of  the  distribution. 

Ahlers^^  compiled  a  data  base  for  a  very  wide  range  of  accidental  and 
test  explosions  and  made  linear  and  quadratic  curve  fits  for  the  data  base, 
including  maximum  debris  distance  versus  yield,  energy  scaled  (W^/^)  maximum 
distance  versus  yield,  and  maximum  debris  distance  versus  impulse.  The 
correlalrion  coefficient  for  the  curve  fits  tended  to  be  low  (0.7  or  less); 
however,  this  may  be  due  to  the  very  wide  range  of  data  used  for  the  fits.  The 
data  base  will  be  discus  ed  later. 

Analytical  methods  for  predicting  debris  range  incorporate  equations  of 
motion  and  i  set  of  assumed  initial  conditions  which  can  include  concrete 
fragment  mass,  launch  angle,  drag  coefficient,  velocity,  lift  coefficient,  and 
attack  angle.  Various  authors  give  relationships  between  debris  range  and 
velocity.  These  studies  were  discussed  earlier.  Those  methodologies  which 
apply  for  R/C  concrete  (i.e.,  make  no  assumptions  such  as  shape  or  density)  are 
listed  in  Table  10. 

Distribution 

Kulesz,  et  al.,^^  give  a  methodology  for  generating  a  hypothetical 
missile  map  using  Figures  2  and  3  and  the  procedure  given  in  Table  11.  This 
procedure  assumes  all  weights  are  distributed  log  normally  in  a  given  interval 
of  range.  Ahlers^^  studied  the  distribution  of  debris  from  the  Pantex 
accident  noted  earlier.  The  reference  included. plots  of  square  feet /concrete 
fragment  and  square  feet/pound  of  debris  versus  range.  This  analysis,  appears  to 
have  assumed  concentric  rings  of  equal  area  about  the  charge  for  the  ground  area 
term.  However,  as  the  reference  points  out,  the  debris  pattern  was 
directional.  Hence,  using  the  area  of  a  whole  ring  instead  of  the  part  of  the 
ring  containing  debris  will  give  lower  ground  area  per  number  or  pound  of 
debris.  This  reference  also  includes  plots  of  square  feet  per  concrete  fragment 
versus  ground  range  for  numerous  accidents  of  a  wide  range  of  structural  types 
(light-frame  construction,  R/C,  and  earth-covered  igloo  magazines).  Merz,  et 
al.,^®  discuss  R/C  debris  distribution  and  the  results  of  1/10  scale  model 
test  for  above-ground  storage  magazines  and  includes  these  conditions: 


FIGURE  3.  CUMULATIVE  PROBABILITY  DISTRIBUTION,  DEBRIS  RANGE  (FT) 


FIGURE  4.  DEBRIU  CUMULATIVE  PROBABILITY  DISTRIBUTION,  NONDIMENS lONAL  RANGE 
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TABLE  10.  METHODOLOGY  FOR  PREDICTION  OF  DEBRIS  RANGE 
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TABLE  10.  (CONT.) 
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24 

TABLE  11.  METHODOLOGY  OF  PREDICTING  DEBRIS  DISTRIBUTION 

A  procedure  for  estimating  the  number  of  debris  missiles  of  a  given 
mass  Interval  which  will  fall  within  a  given  distance  from  an  explosion 
inside  a  building  is  as  follows : 

1.  Estimate  =  total  destroyed  weight  of  the  building  (portion 
of  the  building  which  has  fragmented).  This  estimate  will  de¬ 
pend  mainly  upon  the  amount  of  explosive  stored  or  machined  in 
the  building  at  an.,  given  time  and  the  building  structure  and 
shape. 

2.  Using  the  weight  distribution  in  Figure  2,  obtain  the  average 
weight  of  debris  from  the  explosion,  W^,  by  reading  it  off 
the  appropriate  curve  at  the  50th  percentile.  The  total  num¬ 
ber  of  debris  fragments  from  the  explosion  is  then 


3.  Using  the  range  distribution  in  Figure  3,  take  equal  percen¬ 
tage  increments  (0-10%,  10-20%,  etc.)  or  equal  range  incre¬ 
ments  (0-10  ft,  10-20  ft,  etc.)  and  find  the  number  of  debris 
fragments  Nj^,  in  each  increment.  (If  equal  percentage  incre¬ 
ments  were  taken,  the  number  of  debris  fragments  in  each  incre¬ 
ment  is.  of  course,  the  same.) 
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TABLE  11.  (CONT.) 

Again  using  the  weight  distribution  In  Figure  2.  determine 
the  percintage  of  debris  In  a  particular  weight  interval.  The 
total  numbers  In  each  range  Interval  have  already  been  calcu¬ 
lated  (Step  3).  Thus,  the  number  of  debris  of  a  particular 
weight  in  a  particular  range  Interval  (distance  out  from  the 
source)  can  be  determined.  The  tiiajor  assumption  made  in  this 
procedure  is  that  all  weights  are  distributed  log  normally  in 
a  given  interval  of  range.  Since  we  could  find  no  correlation 
between  weight  and  range  for  a  given  energy  level,  and  since 
weight  is  log  normally  distributed  over  each  energy  level  (which 
covers  the  entire  range) ,  there  is  no  reason  to  assume  that 
weight  is  not  log  normally  distributed  within  a  given  range 
Increment. 
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o  Debris  dispersion  in  R/C  buildings  is  dependent  on  building  shape. 
In  tb-*  direction  normal  to  walls  there  is  a  greater  number  of 
concrete  fragments  with  larger  range  than  out  from  corners. 

o  iTie  more  strength  and  ductility  in  a  structure,  the  more 

d ■■  rectional  the  dispersal  is  in  the  surface  normal  directions.  The 
r.oc'il  magazines  were  reinforced  concrete  which  represented  ductile 
consiruction  of  considerable  strength  and  the  failure  of  the 
maga:’’ nea  were  similar  to  a  chamber  pressure  failure,  with  forces 
actl.ig  perpendicularly  on  the  four  walls.  The  initial  shock  waves 
did  not  contain  enough  energy  to  fail  the  structure. 

b  Merz  saggasts  that,  because  of  the  relatively  flat  trajectories  of 
fragments,  distributions  should  not  be  measured  by  number  of 
concrete  fragments  per  ground  area,  but  number  of  concrete 
fragments  transversing  a  plane  vertical  to  the  ground  and  parallel 
to  a  surface.  This  is  because  debris  presents  a  hazard  along  its 
entire  path  and  not  just  where  it  eventually  lands. 

o  Merz  also  pioints  out  that  for  the  model  tests  and  many  full-scale 
instances,  the  debris  presents  the  predominate  hazard  in  an 
accident  of  this  type. 


Final  velocity  of  debris  fragments  can  be  predicted  using  equations  of 
motion,  the  initial  conditions  and  solving  for  conditions  at  impact.  These 
equations  are  given  in  Table  10. 


Hazards  Analysis 


Moseley  and  Whitney^^  combined  methodologies  in  several  of  the  above 
discussed  categories  for  a  prediction  of  a  debris  hazard  for  an  explosion  in  an 
R/C  structure.  The  object  of  the  project  was  to  investigate  several  analytical 
methods  of  predicting  blast  and  fragment  hazards  from  an  accidental  explosion 
inside  a  Norwegian  aircraft  shelter.  This  work  was  in  support  of  model  tests 
being  conducted  in  Norway.  The  steps  taken  for  the  fragment  analysis  included: 


o  Determination  of  debris  fragment  size  and  initial  trajectory 

angle.  Films  of  model  tests  were  reviewed  to  obtain  predominate 
launch  angles  and  relationship  between  debris  fragment  sizes  and 
location  on  structure.  Trajectories  were  found  to  be  approximately 
normal  to  the  c. rface.  Debris  sizes  were  found  to  be  small,  close 
to  the  off-centered  (in  the  building)  charge  and  a  wide  range  of 
sizes,  including  entire  surfaces,  away  from  the  charge  location. 

o  Determination  of  initial  velocity  consisted  of  calculating  velocity 
from  both  shock  loading  and  quasi-static  pressures.  Shock 
velocity,  V^,  was  calculated  using  Equation  (4).  Quasi-static 
velocity  was  calculated  using  the  program  CYLIN.  The  velocities 
computed  using  both  shock  and  quasi-static  loads  appeared  to 
compare  well  with  the  model  tests.  The  structure  was  overloaded 
severely  in  these  tests. 


37 


NSWC  TR  85-114 


o'  Debris  range  was  determined  using  the  program  FRISB,  Input  was 
concrete  fragment  sizes  and  mass,  and  corresponding  velocities  as 
determined  above.  Launch  angles  were  predominately  normal  to  the 
surfaces,  but  also  included  some  off-normal  launch  angles  indicated 
by  films.  The  code  was  run  considering  only  drag  forces,  then 
repeated  considering  both  lift  and  drag  forces.  These  calculations 
led  to  determination  of  maximum  concrete  fragment  ranges  as  a 
function  of  direction  around  the  shelter.  Comparison  with  the 
results  of  the  test  data  are  incomplete  pending  further  test  and 
data  reduction. 

Data  Bases 

Several  reports  include  the  bulk  data  or  indicate  where  bulk  data  for  R/C 
debris  may  exist.  The  data  bases  available  include: 

o  Kulesz,  et  al.,^^  include  debris  data  from  six  structures 

obtained  from  DDESB  accident  reports.  The  data  include  estimated 
yield,  individual  debris  weight,  and  individual  debris  range. 

o  Reference  25  contains  a  large  collection  of  accident  data  which 
include  amount  of  explosive,  maximum  debris  distance,  assumed 
applied  impulse,  and  W^' ^-scaled  debris  distance  and  impulse. 

Also  the  reference  includes  a  table  of  data  from  the  Pantex 
accident  including  individual  debris  size,  mass,  and  range.  Number 
of  concrete  fragments,  total  mass  of  debris,  average  debris  weight, 
and  individual  debris  mass  are  given  for  ground  range  increments  in 
this  reference. 

o  References  27  and  28  indicate  data  bases  for  modal  tests  are 
available,  however,  do  not  include  the  data  in  these  reports. 

o  Reference  26  includes  the  data  base  for  debris  from  scale  models  of 
dividing  walls.  This  includes  number  of  concrete  fragments 
recovered  for  several  mass  ranges,  average  mass  for  each  range, 
average  distance  for  each  mass  range,  total  number  of  concrete 
fragments,  and  concrete  fragment  shape.  The  reference  includes  a 
model  analysis  and  plots  of  scaled  data. 

o  Edmunds^^  has  debris  data  for  brick  walls  which  fail  under 

instantaneously  applied  static  loads.  This  includes  quantity, 
weight,  and  distance  of  individual  debris.. 
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CHAPTER  3 


ACCEPTABLE  HAZARD  HANDLING  ARC  CRITERIA 


Currently,  manufacturing,  storage,  or  handling  of  high  explosives  in 
quantities  between  68  and  13,600  kg  (150  and  30,000  lb)  is  regulated  by  the 
Explosive  Safety  Quantity  Distance  Criteria  established  by  the  DDESB.  These 
criteria  define  the  Acceptable  Hazard  Handling  Arc  as  381  m  (1250  ft)> 

Exemptions  may  be  granted  for  those  cases  where  approved  analysis  or  testing  has 
been  conducted.  In  this  case,  the  acceptable  hazard  handling  arc  for  an 
explosion  event  is  defined  by  the  minimum  range  at  which  both  the  blast 
overpressure  and  the  fragment  hazard  criteria  are  satisfied.  The  criteria  are 
as  follows: 

1.  the  blast  overpressure  is  less  than  6.89  kPa  (1.0  psi),  and 

2.  the  number  of  hazardous  fragments  entering  a  given  region  is  less 

than  1.0  hazardous  fragment  per  55.7  (600  ft^)  of  ground 

surface  area.  A  fragment  is  considered  hazardous  when  it  has  an 
impact  energy  equal  to  or  in  excess  of  8UJ  (58  ft-lb). 

During  this  program  we  attempted  to  discover  the  basis  for  the  above  criteria. 

We  did  not  identify  the  original  source,  but  we  believe  the  8UJ  (58  ft-lb) 
refers  to  a  50  percent  injury  (50  percent  fatality)  level  of  a  combat  soldier, 
and  the  55.7  m^  (600  ft^)  refers  to  100  times  the  ground  surface  covered  by 
one  prone  soldier.  In  Reference  30,  a  review  of  different  criteria  for 
hazardous  fra^ents  was  presented.  These  different  criteria  are  summarized  in 
Table  12.^^*^*  Two  of  the  criteria  are  based  on  the  kinetic  energy  at  impact, 
while  the  other  criterion  is  based  on  impact  momentum.  Interestingly  enough, 
the  last  three  criteria  are  all  less  conservative  than  the  standard  80J  (50 
ft-lb)  criteria.  This  is  shown  in  Table  13,  where  allowable  velocities  for 
fragments  with  masses  of  0.025,  O.l  and  0.4  kg  (0.055,  0.22,  and  0.90  lb)  are 
summarized  for  all  four  criteria.  If  the  last  three  criteria  were  used  instead 
of  the  80J  (58  ft-lb)'  criteria,  then  larger  acceptable  hazard  arcs  would  be 
derived. 


The  reason  for  this  discussion  is  to  emphasize  that  different  definitions 
for  a  hazardous  fragment  exist  in  the  literature,  and  for  some  situations  a 
different  criterion  than  80J  (58  ft-lb)  of  energy  may  be  more  appropriate. 
Indeed  a  more  rational  approach  to  specifying  hazard  arcs  would  be  to  specify 
the  fragment  energy  as  a  function  of  the  "target.”  As  an  example,  different 
hazardous  fragment  energies  could  be  derived  for  a  specified  level  of: 
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TABLE 

12.  HAZARDOUS  FRAGMENT  CRITERIA 

Criterion  , 

Definition 

Reference 

Casualty  Level 

A 

K.E. 

=  80J 

Wide  use 

B 

K.E. 

D 

=»  2.5  X  lO^N 

Netherlands^® 

Penetration  of 
37  mm  of  poplar 

C 

=  27.1  kg^'^m/s 

German,  1965^^ 

D 

MV 

A 

-  2.0  X  10^  ^ 

m-s 

32 

U.S.,  1951 

K.E.  -  fragment  kinetic  energy 

M  -  fragment  mass 

V  fragment  velocity 

A  -  average  fragment  presented  area 

D  -  characteristic  fragment  length 
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TABLE  13.  ALLOWABLE  VELOCITIES  FOR  CERTAIN  MASSES  USING  THE 
HAZARDOUS  FRAGMENT  CRITERIA 


Criterion 

Allowable 
0.025  kg 

Velocities  (m/s) 
0.1  kg 

for  Fxagment  Masses  of 
0.4  kg 

A 

80 

40 

20 

B* 

190 

120 

75.8 

C 

119 

68.1 

39.1 

* 

D 

204 

132 

83.1 

Assuming  spherical,  steel  fragments 


41 


NSWC  TR  83-114 


1.  human  casualty,! 

2.  damage  to  nearby  buildings  (concrete  structures),  and 

3.  damage  to  nearby  ships  (steel  structures). 

However,  for  safety  evaluations,  the  human  casualty  level  would  probably  always 
override  the  damage  levels  for  buildings  and  ships.  Correspondingly  the  area 
used  in  the  estimation  of  the  fragment  flux  should  also  correspond  to  the 
appropriate  presented  area  for  a  human,  a  building,  or  a  ship.  In  the 
estimation  of  areas,  the  presented  area  should  correspond  to  the  frontal  area 
for  low  trajectory  fragments,  or  ground  surface  area  for  high  trajectory 
fragments.  For  safety  evaluations,  the  soldiers  position  should  be  standing 
(not  prone  as  is  now  considered) — for  this  case,  then,  the  vulnerable  target 
area  for  low  trajectory  debris  is  greater  than  for  high  trajectory  debris.  The 
use  of  more  specific  definitions  for  hazardous  fragment  flux  would  result  in 
more  credible  safe  handling  arcs. 


A  GENERAL  METHOD  FOR  PREDICTING  IHE  FRAGMENT  AND  DEBRIS  ENVIRONMENT 

The  general  method  for  predicting  the  fragment  and  debris  environment 
resulting  from  an  accidental  explosion  is  given  in  flow  chart  form  in  Figure  5. 
Tile  method  is  general  and  can  be  applied  to  an  explosion  of  a  weapon  in  the  open 
field  or  to  an  explosion  inside  a  building.  For  the  case  of  an  explosion  in  a 
building,  the  methodology  is  applied  separately  for  the  fragments  and  the 
debris,  and  the  results  are  combined  in  the  last  step  to  define  the  hazard  arc. 
The  method  can  conceptually  be  applied  to  situations  where  no  data  or  limited 
data  are  available  by  substitution  of  engineering  estimates  for  the  various 
parameters  in  the  problem  definition  phase. 

The  methodology  presented  in  this  interim  report  parallels  a  similar 
methodology  developed  by  J.  M.  Ward^^  for  predicting  and  subsequently 
analyzing  the  fragment  hazard  associated  with  the  accidental  detonation  of 
weapons,  specifically  a  pallet  load  of  Mk  82  bombs.  His  methodology,  which  is 
an  application  of  concepts  developed  by  F.  B.  Porzel^^  for  NESIP,  was 
developed  for  the  special  case  of  an  explosion  in  the  open  field.  This  method 
provides  for  both  high-  and  low-launch  angle  fragments.  This  distinction  is 
important  in  the  establishment  of  acceptable  hazard  arcs.  In  this  report,  we 
have  included  this  important  concept  in  the  estimation  of  fragment  hazards  for 
both  fragments  and  debris  missiles. 

Assumptions 

Several  assumptions  were  made  in  the  preparation  of  this  flow  chart.  These 
assumptions  are: 

1.  There  is  a  uniform  distribution  of  fragments  or  debris  missiles 
with  respect  to  launch  angle. 

2.  At  any  launch  angle,  there  is  a  uniform  distribution  of  fragment  or 
debris  mass. 
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FIGURE  5.  FLOW  CHART  FOR  PREDICTING  THE  FRAGMENT  AND  DEBRIS 
ENVIRONMENT  AROUND  EXPLOSIONS 


/ 
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3,  The  fragments  or  debris  have  the  same  initial  velocity. 

4.  The  fragments  or  debris  missiles  crhibit  drag  and  no  lift.  The 
drag  area  is  approximated  from  a  j^mple  function  of  the  fragment 
mass.  The  function  is  derived  from  measurements  of  the  fragment 
presented  area  for  a  small  sample  of  fragments.  The  drag 
coefficient  is  chosen  according  to  the  fragment  shape.  Drag 
coefficients  for  tumbling  cubes  or  for  shell  fragments  are  used  in 
the  absence  of  specific  fragment  shape  informat  on. 

The  first  three  assumptions,  which  are  related,  are  made  to  simplify  the  problem 
formulation.  It  is  realized  that  in  practice  uniform  distributions  are  not 
likely  to  be  encountered.  However,  the  inclusion  of  nonuniform  distributions 
(where  they  can  be  derived)  would  require  Monte-Carlo  simulations  in  order  to 
arrive  at  a  reasonable  solution.  For  some  situations,  such  as  the  off-center 
explosion  of  a  charge  within  a  building,  the  assumption  of  a  constant  initial 
velocity  may  not  be  appropriate .  In  this  case,  it  is  possible  to  use  this 
methodology  individually  for,  as  an  example,  each  of  the  four  walls  and  the 
roof,  and  to  superimpose  the  results  in  the  last  stage  of  the  method.  The  last 
assumption  refers  I  Co  the  method  of  calcul  -  ion  used  to  estimate  the  debfis 
trajectory.  Drag  coefficients  and  areas  for  stable  orientation  solids  are  used 
in  the  trajectory  calculations,  since  little  information  is  available  for 
equivalent  drag  coefficients  for  tumbling  objects. 

APPLICATION  OF  THE  METHOD 
Problem  Definition 

The  fragment/debris  environment  prediction  methodology  presented  in 
Figure  5  consists  of  four  parts.  The  first  and  moit  difficult  part  to  calculate 
consists  of  the  problem  definition.  Where  data  are  available,  this  phase, 
although  tedious,  presents  little  problem.  Where  no  data  are  available,  a 
considerable  amount  of  engineering  judgement  is  required  to  define  adequately 
and  accurately  the  requisite  input  parameters.  The  specific  parameters  which 
must  be  defined  in  Che  initial  part  of  the  methodology  are: 

1.  Che  mass  distribution — generally  presented  in  the  Mott  format: 

N(m  >  M)  =  Nq  exp  (-m/nio)^ 

2.  Che  initial  velocity  of  the  fragments  or  the  debris 

3.  the  range  of  launch  angles — for  example  debris  missiles  from  a 
reinforced  concrete  wall  preferentially  emanate  normal  to  the 
original  plane  of  the  wall.  Thus  the  range  of  launch  angles  can  be 
fairly  narrow. 

4.  the  fragment/debris  drag  areas  aind  coefficients — generally  tumbling 
cubes  and  shell  shapes  are  selected  for  the  fragments,  and  drag 
coefficients  are  taken  from  standard  tables.  The  drag  areas  are 
estimated  from  a  simple  relationship  of  presented  area  as  a 
(unction  of  the  fragment  mass.  For  debris,  other  shapes  may  be 
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more  appropriate;  for  example,  concrete  rubble  might  be  represented 
by  tumbling  cubes  or  spheres. 

5.  the  range  distribution—the  number  of  fragments/debris  per  recovery 
zone  broken  down  by  mass,  N  "  f(R,m),  determined  experimentally.* 

In  this  section  tite  methodology  will  be  illustrated  using  data  collected  in 
some  recent  experiments  conducted  at  SwRI.26  in  these  tests,  reinforced 
concrete  panels  were  overloaded  by  bare  explosive  charges  placed  at  small  scaled 
distances  from  the  panel.  In  each  test,  the  debris  was  recovered,  and  the 
recovery  location,  mass,  shape,  and  physical  dimensions  were  recorded.  The 
test,  presented  as  an  example  in  this  report,  is  summarized  in  Table  14.  A 
missile  map  for  this  test  is  given  in  Figure  6,  and  a  mass  distribution  is  given 
in  Figure  7.'  The  maximum  velocity  measured  in  the  high-speed  films  was  23.5  m/s 
(77  fps).  The  number  of  concrete  fragments  collected  in  the  six  recovery  areas 
marked  on  the  missile  map  is  summarized  in  Table  15. 

The  range  of  launch  angles  assumed  by  the  debris  missiles  was  established 
based  on  the  missile  map  given  in  Figure  6.  Assuming  that  the  debris  is  ejected 
in  a  circularly  symmetric  pattern,  the  range  of  launch  angles  is  the  same  as  the 
range  of  the  debris  spray  angle.  This  assumption  was  verified,  at  least 
approximately  by  the  high-speed  films  taken  during  the  test.  It  was  noted  in 

the  films  that  the  preferential  debris  launch  angle  in  this  test  was  normal  to 

the  original  plane  of  the  wall  panel,  which  is  in  conflict  with  the  assumption 
of  a  uniform  distribution  of  debris  missiles  with  respect  to  launch  angle. 
Nonetheless,  for  the  purposes  of  illustrating  the  methodology,  the  range  of 
launch  angles  used  in  this  example  problem  was  taken  to  be  -20  to  +80  degrees. 
Each  debris  missile  was  assumed  to  have  an  initial  height  of  9  inches,  the 
height  of  the  charge  above  the  ground.  This  represents  another  simplification 
for  the  purpose  of  this  example  problem.  For  the  test,  of  course,  the  initial 

height  varied  from  ground  level  to  the  top  of  the  panel. 

The  final  input,  which  was  necessary  to  demonstrate  the  methodology,  was 
debris  drag  or  presented  area.  To  generalize  the  effects  of  different  debris 
shapes  (pancake  or  chunky  type)  and  mode  of  flight  for  tumbling  concrete 
fragments,  average  areas  were  determined  for  each  debris  mass  recovered  by 
computing  the  equivalent  volume  of  the  debris,  using  the  three  orthogonal 
dimensions  corresponding  somewhat  to  width,  length,  and  thickness,  and  assuming 
the  debris  missile  was  a  rectangular  solid.  This  leads  to  the  following 
expression  for  the  drag  area; 

A  =  (L  X  W  X  T)2/3  (5) 

The  areas  calculated  using  Equation  (5)  were  then  plotted  as  a  function  of 
debris  mass  as  shown  in  Figure  8.  The  equation  of  this  line,  the  percent 
standard  deviation  and  the  multiple  correlation  coefficient,  R,  are  included  on 


*This  parameter  is  considered  as  an  input  variable  only  -if  there  are  data 
■available.  When  test  data  are  available  then  the  analytical  procedure  described 
is  a  procedure  for  analyzing  the  recovered  f ragment/debris  data. 
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FIGURE  6.  MISSILE  MAP  FOR  A  TYPICAL  DIVIDING  WALL  TEST 
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TABLE  15.  NUMBER  OF  CONCRETE  FRAGMENTS  RECOVERED  BY  MASS 
IN  THE  SIX  RECOVERY  ZONES 

Number  of  Debris  Missiles  Recovered  in  Zone 
Mass  Interval  (gm) _ 1 _ 2 _ 3 _ 4  5 _ 6 

0.0  -  0.9  4  0  0  1  1  0 

0.9  -  4.5  20  3  0  17  8  0 


4.5  -  27.0  13  94  4,4  3 

27.0  -  90.0  100  100 


Che  graph  in  Figure  8.  Using  this  fitted  equation,  drag  areas  could  be 
determined  for  all  masses  considered  in  Che  analysis.  An  average  drag 
coefficient  of  0.9  was  used  for  all  debris  in  this  analysis  since  they  were 
mostly  chunky  in  shape  and  0.9  is  the  average  between  the  drag  coefficient  for  a 
face-on  and  edge-on  cube. 

Hazard  Arc  Definition 

General.  In  this  report,  two  criteria  for  defining  debris  as  hazardous  to 
personnel  are  discussed.  The  traditional  criterion  considers  debris  hazardous 
if  it  strikes  Che  ground  surface  with  more  Chan  80J  (58  fc-lb)  of  kinetic 
energy.  This  criterion  applies  to  prone  personnel.  Another  criterion  considers 
debris  hazardous  at  a  specific  point  if  it  passes  through  a  vertical  plane  with 
more  than  80J  (58  ft-lb).  This  criterion  applies  to  standing  personnel  and, 
therefore,  the  height  of  the  vertical  plane  is  set  at  3  m  (9.8  ft).*  The 
following  subsection  will  describe  the  methodology  which  will  be  used  to 
calculate  debris  densities  in  the  horizontal  and  vertical  planes  around  an 
explosive  source.  This  methodology  requires  that  Che  ground  surface  around  Che 
explosion  source  be  divided  into  recovery  zones.  The  horizontal  debris  density 
in  each  recovery  area  is  calculated  as  Che  number  of  hazardous  debris  elements 
Chat  land  in  that  zone  divided  by  its  ground  surface  area.  The  vertical  debris 
density  corresponding  Co  Chat  recovery  area  is  calculated  as  the  number  of 
hazardous  debris  elements  that  cross  any  vertical  plane  3m  (9.8  ft)  high  in  the 
recovery  zone.  One  such  vertical  plane  is  shown  schematically  in  Figure  9.  A 
recovery  area  is  considered  unsafe  for  personnel  if  the  hazardous  debris  density 
in  either  the  vertical  or  horizontal  plane  exceeds  one  hazardous  debris  element 
per  55.7  m^  (1/600  ft^). 

Debris  Density  in  the  Horizontal  Plane.  The  hazardous  debris  in  the 
horizontal  plane  is  determined  by  examining  the  debris  which  land  in  each 
recovery  area  in  a  test.  For  each  recovery  location  there  are  two  trajectory 
solutions  which  could  explain  the  observed  debris  range:  a  high  and  a  low 
trajectory.  This  distinction  is  important  since  the  input  energies  associated 
with  these  two  trajectory  solutions  can  be  remarkably  different.  For  this 
reason  it  is  necessary  to  determine  the  ratio  of  high  trajectories  to  low 
trajectories  per  recovery  area.  To  accomplish  this,  typical  or  representative 


*The  height  of  an  average  man  plus  a  50  percent  safety  margin. 
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FIGURE  9.  DEFINITION  OF  VERTICAL  PLANE  USED  TO  ESTIMATE 
THE  VERTICAL  DEBRIS  DENSITY 
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masses  and  launch  angles  are  selected  and  the  debris  range  is  calculated  using  a 
standard  trajectory  code.*  In  this  illustrative  example,  the  calculations  were 
performed  for: 

masses:  0.9,  4.5,  27.0,  90  gm  (0.002,  O.OI,  0.06,  0.20  lb) 

launch  angles:  -20  to  80  degrees  in  10  degree  increments 
initial  velocity:  23.5  m/s  (77.1  fps) 

drag  area:  determined  from  Figure  8  and  the  debris  mass 

By  systematically  varying  the  launch  angles,  it  is  possible  to  define  the 
maximum  possible  range  as  a  function  of  the  launch  angle  for  each  mass  group. 

One  such  plot  is  given  in  Figure  10.  This  figure  is  used  to  define  the  low 
trajectories,  in  this  case  all  trajectories,  resulting  from  a  launch  angle  less 
than  42  degrees.  Similarly  high  trajectories  are  defined  as  having  launch 
angles  greater  than  42  degrees. 

Figure  10  can  be  used  to  establish  the  proportion  of  low  to  high  trajectory 
debris  with  a  mass  of  90  gm  (0.20  lb)  recovered  in  a  particular  zone.  Consider 
recovery  zone  5  which  is  12  to  18  m  (20  to  40  ft)  downrange.  Debris  with  a  mass 
of  90  gm  (0.20  lb)  could  have  landed  in  this  recovery  zone  if  it  was  launched  at 
an  initial  angle  between  6.7  and  10  degrees  or  between  76.7  and  81.7  degrees. 

The  range  of  possible  low  trajectory  launch  angles  is  called  and  is 
10  -  6.7  or  3.3  degrees  wide.  The  range  of  high  trajectory  angles  is  ASy, 
which  f'-t  this  case  is  5.0  degrees  wide.  The  percentage  of  low  trajectory 
debris  elements  in  the  recovery  area  is  then  defined  as: 


49, 


L  *  AOjj 


0.398 


(6) 


The  percentage  of  high  trajectory  debris  elements  is  then: 


Fh  -  1.0  -  Fl  -  0.602 


(7) 


Once  the  proportion  of  high  and  low  angle  trajectories  are  known  for  each 
mass  group,  the  number  of  hazardous  debris  in  the  recovery  zone  can  be 
calculated.  This  is  accomplished  by  examining  the  trajectory  code  output  for 
each  representative  mass  and  an  average  high  trajectory  launch  angle, 

If  the  kinetic  energy  of  this  case  exceeds  the  80J  (58  ft-lb)  criterion,  then 
the  number  of  hazardous  debris  elements  in  that  recovery  area  is  incremented  by: 

%  »  Nh  +  N(m)  •  fji(m)  (8) 


*At  SwRI,  the  trajectory  code  used  is  FRISB,  which  was  developed  under  a  NASA 
contract. 
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FIGURE  10.  MAXIMUM  RANGE  ENVELOPE  FOR  90  GM  (0.20  LB)  DEBRIS 
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where 

%  ■  number  of  hazardous  debris  elements 

N(m)  v  number  of  debris  elements  with  a  mass  m  found  in  the  recovery 
zone 

fl}(m)  ■  percentage  high  trajectory  debris  elements  with  a  mass 
m  in  the  recovery  zone 


Similarly,  the  contribution  of  low  angle  trajectory  debris  of  the  specific 
mass  group  to  the  total  number  of  hazardous  debris  in  the  recovery  can  be 
estimated  by  examining  the  trajectory  code  output  again  for  each  representative 
mass  and  an  average  low  trajectory  launch  angle,  S  Whenever  the  impact 
kinetic  energy  of  one  of  these  cases  exceeds  the  kinetic  energy  criterion,  then 
the  number  of  hazardous  debris  in  that  recovery  area  is  incremented  by: 

%  “  Nh  +  N(m)  •  f^Cm)  (9) 

The  methodology  for  the  calculation  of  the  hazardous  debris  density  was 
applied  to  the  dividing  wall  test  data.  The  results  of  these  calculations  are 
summarized  in  Table  16.  For  each  of  the  six  recovery  zones,  the  four 
representative  masses  and  the  high  and  low  angles  (9}^  and  9j{),  the  code 
was  examined  and  the  debris  kinetic  energy  at  impact  was  calculated.  Using  a 
properly  scaled  hazard  criterion  of  0.37J  (0.27  ft-lb),*  the  number  of  hazardous 
debris  elements  was  determined  for  each  recovery  zone.  The  hazardous  debris 
density  was  calculated  and  normalized  to  55.7  (600  ft^).  The  results 

summarized  in  Table  16  indicate  that  the  minimum  hazard  arc  is  beyond  recovery 
zones  3  and  6  or  more  than  24  m  (78  ft)  from  ground  zero  (see  Figure  6). 

It  should  be  noted  that  for  some  classes  of  events,  the  hazard  criterion 
(horizontal  plane)  will  be  satisfied  at  points  close  to  ground  zero,  but  will 
not  be  satisfied  at  large  distance  from  ground  zero.  This  result  can  be  obtained 

*The  dividing  wall  test  program  was  conducted  at  1/6  scale.  Well  established 
laws  of  scaling  require  that  energies  scale  as  the  geometric  scale  factor 
cubed.  Thus  the  energy  criterion, of  80J  (58  ft-lb)  becomes  0.37J  (0.27  ft-lb). 
The  range  that  a  debris  raissilt  "travels  is  not  properly  scaled  in  the  model 
tests  since  the  gravity  field  w  s  not  increased  by  the  geometric  scale  factor 
(1/6)  in  the  tests.  This  means  the  model  scale  range  observations  are  of  the' 
same  magnitude  as  those  observed  in  full  scale.  Therefore,  the  debris  density 
of  1/55.7  m^  (1/600  ft^)  is  appropriate.  Note  that,  neglecting  drag,  the 
range  of  a  debris  element  is; 


R 


V  ^  sin  29 
o 


g 


+  V 
o 


d  is  the  initial  height  of  the  debris. 


In  these  experiments,  the  velocity,  Vq,  and  the  gravity,  g,  were  held 
constant  and  the  initial  height  of  the  debris  was  properly  scaled.  If  we 
consider  debris  emerging  from  the  ground  surface  (d  *  O),  the  model  and 
full-scale  ranges  are  identical.  However,  for  proper  scaling,  the  model  ranges 
should  be  modified  to  account  for  the  scale  change  of  1/6  for  the  second  term 
in  this  range  equation. 
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TABLE  16.  CALCULATIONS  OF  THE  HAZARDOUS  I 


Recovery  Zone 


*1.  »>■ 


NaAfi 

(«■) 


1 


3.7 


♦  .5 
27.0 
90,0 


2 


8.S 


0.398  0.9 

4.S 
27.0 
90.0 


79.2  0.602  0.9 

6.5 
27.0 
90.0 


3 


13.2  0.562  0.9 

6.5 
27.0 
90.0 

76.0  0.658  0.9 

4.5 
27.0 
90.0 


DENSITY  IN  THE  HORIZONTAL  PLANE 


X 

KFf 

(J) 

M(») 

"h 

15.2 

C.105 

6 

- 

18.0 

0.736 

20 

20 

20.2 

5.6 

13 

13 

21.3 

20.7 

1 

1 

Tplal  M 

In  Zone  1 

34 

Razarjous  Density  * 

206 

12.5 

0.071 

0 

- 

15.5 

0.566 

3 

1.19 

18.4 

6.62 

9 

3.58 

19.9 

18.0 

0 

- 

8.9 

0.036 

0 

- 

11.6 

0.292 

3 

- 

14.7 

2.92 

9 

5.42 

16.7 

12.6 

0 

- 

Total  Ml, 

in  Zone  2 

10.2 

Harardo->s 

Density 

61.1 

10.9 

0.056 

0 

- 

13.8 

0.636 

0 

- 

16.9 

3.91 

4 

2.17 

18.7 

15.9 

0 

- 

10.5 

0.050 

0 

- 

13.5 

0.610 

0 

- 

16.6 

3.72 

4 

183 

18.3 

15.1 

0 

- 

Total  Nji 

in  Zone  3 

4.00 

Hazardous 

Density 

26.0 
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TABLE  16.  (CONT.) 


(deg) 

Haas 

Vf 

KEf 

N(i«) 

Recovery  Zona 

(go) 

(■/a) 

(J) 

4 

3.7  1.0 

0.9 

15.2 

0.105 

1 

- 

4.5 

18.0 

0.  734 

17 

17 

2T.0 

20.2 

5.6 

4 

'  ■  4 

90.0 

21.3 

20.7 

1 

_1 _ 

1 

Total 

in  Zone  4 

22 

Hazardous 

Density 

132 

5 

8.5  0.398 

0^.9 

12.5 

0.071 

1 

- 

4.5 

15.5 

0.546 

8 

3.18 

27.0 

18.4 

4.62 

4 

1.59 

90.0 

19.9 

18.0 

0 

“ 

79.2  0.602 

0.9 

8.9 

0.036 

1 

. 

4.S 

11.4 

0.292 

8 

- 

27.0 

14.7 

2.92 

4 

2.41 

90.0 

16.7 

12.6 

0 

- 

Total  Njj 

In  Zone  5 

7.18 

Hazardous 

Density 

43.1 

6 

13.2  0.542 

0.9 

10.9 

0.054 

0 

- 

4  ..5 

13.8 

0.434 

0 

- 

27.0 

16.9 

3.91 

3 

1.63 

90.0 

18.7 

15.9 

0 

- 

0.458 

0.9 

10.5 

0.050 

0 

- 

4.5 

13.5 

0.410  ' 

0 

,  - 

27.0 

16.6 

3.72 

3 

1.37 

90.0 

18.3 

15.1 

0 

- 

Total  Nji 

in  Zone  6 

3.00 

Hazardous 

Density 

18.0 

*  Number  of 

hazardous  dehrla  divided  by  ' 

9.29  (100  ft^) 

ground  surf 

ace  area 

and  normnll^-cd  to  55.7  n  (600  ft  ) 

t  In  tills  Illustrative  example,  the  f  and  f  used 

correspond  to  the 

90 

(0.20  lb) 

deb^^*;  rlement.  In.  practice 

,  separate 

estimates  for  f  and 

should  he 

mad^  for  each  mass  and  range 

group. 
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when  the  debris  missiles  emerge  predominately  with  launch  angles  close  to  the 
optimum  value  for  maximum  range.  Thus  the  safety  analyst  must  be  careful  to 
analyze  the  full  range  of  plausible  launch  angles. 

Debris  Density  in  the  Vertical  Plane.  The  debris  density  calculations  for 
the  vertical  plane  parallel  those  for  the  horizontal  plane.  The  difference  is 
that  the  debris  velocity  and  its  vortical  height  as  it  crosses  a  vertical  plane 
are  used  (instead  of  the  range  and  the  velocity  at  impact)  in  the  determination 
of  whether  it  is  hazardous.  In  practice,  all  combinations  of  the  representative 
debris  masses  and  the  representative  launch  angles  are  used  in  the  trajectory 
calculations.  As  each  debris  element  crosses  the  recovery  area,  its  velocity 
and  height  are  noted.  If  the  debris  height  is  less  than  3  m  (9.8  ft),  and  the 
debris  kinetic  energy  exceeds  80J  (58  ft-lb)  then  that  piece  of  debris  is 
considered  hazardous  and  the  number  of  hazardous  debris  in  that  recovery  zone  is 
incremented  according  to: 


Nv  “  Ny  +  Fi^(m)N(m)  (10) 

where 

My  «  number  of  hazardous  debris  in  the  vertical  plane 

FL(m)  *  percentage  low  trajectories  for  a  mass  group 
(see  Equation  (6)) 

N(m)  =  number  of  debris  elements  with  mass,  m,  which  landed 
in  the  recovery  zone  being  considered 

When  all  of  the  debris  has  been  accounted  for,  the  hazardous  debris  density  in 
the  vertical  plane  is  estimated  by  dividing  My  by  the  area  of  the  vertical 
plane  at  the  midspan  of  the  recovery  area.  A  recovery  area  is  considered  unsafe 
for  personnel  if  the  hazardous  deb -is  areal  density  exceeds  1/55.7  m^,  (1/600 
lt2>. 


As  an  illustrative  example,  consider  again  the  dividing  wall  test  data. 
Trajectory  calculations  were  performed  for  the  same  masses  and  launch  angles 
considered  for  the  horizontal  debris  density  estimates.  In  this  case  the 
trajectory  code  was  modified  to  print  out  the  debris  velocity  and  height  as  it 
crofsed  vertical  planes  at  the  midpoint  of  the, recovery  areas  located  at  9,  15, 
and  21  m  (30,  50,  and  70  ft)  downrartge  (see  Figure  6).  These  calculations  are 
summarized  in  Table  17.  Using  the  scaled  hazardous  kinetic  energy  of  0.37J 
(0.27  £t-lb),  the  number  of  hazardous  debris  missiles  per  recovery  plane  wis 
calculatei  and  the  hazardous  debris  densities  were  compared  to  1/55.7  m^ 

(1/600  ft2).  The  comparison  shown  at  the  bottom  of  Table  17  indicates  that 
the  minimum  hazard  arc  from  the  vertical  plane  calculations  is  less  than  21  m 
(70  ;c).  Since  this  arc  is  less  than  that  for  the  horizontal  plane 
calcul itions ,  the  overall  safety  arc  is  controlled  by  the  horizontal  plane 
calculations.  It  should  be  noted  that  these  findings  resulted  because  the 
vertical  ylane  calculations  considered  the  plane  located  at  the  midspan  of  each 
of  tie  six  recovery  zones,  not  at  the  close-in  boundary  of  the  recovery  area 
sector.  In  practice,  the  hazard  arc  based  on  the  vertical  plane  calculations 
will  control  the  definition  of  the  minimum  acceptable  hazard  arc. 
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TABLE  n.  NUMBER  OF  HAZARHOUS  DEBRIS  IN  THE  VERTICAL  PLANE 


• 

(det) 

■mat 

(»i) 

Velocity 

9m 

(W.)  at  ch.  Cbrl.  Acia..  ?la.  at  C2) 
13.  21.  1 

(a) 

at  Each  Recovery  Plane 

9a  '  15«  21a 

3.7 

1.0 

0.9 

-<3) 

. 

. 

4 

• 

- 

4.3 

- 

- 

- 

20 

- 

- 

- 

(Zone  1) 

27.0 

20.2 

- 

- 

13 

13 

- 

90.0  1 

21.2 

- 

1 

1 

- 

.  - 

8.3 

0.398 

0.9 

- 

- 

mm 

t 

- 

- 

4.3 

17.3 

- 

- 

1.19 

- 

- 

(Zona  2) 

27.0 

- 

'  18.3 

- 

1 

3.33 

- 

90.0 

• 

20.0 

- 

H 

- 

* 

• 

13.2 

0.342 

0.9 

- 

- 

- 

- 

- 

4.3 

- 

U.3 

- 

0 

- 

- 

- 

(ZoM  3) 

27.0 

- 

- 

- 

4 

- 

- 

- 

90.0 

- 

- 

18.7 

0 

• 

- 

- 

Zoos  1 

Zone  2 

Zoo.  3 

Humber  of  Raaardoua  Debris 

13.19 

3.33 

0 

Dsbrla  Density  '  * 

182.0 

42.4 

0 

3.7 

1.0 

0.9 

mm 

- 

• 

1 

. 

- 

* 

4.5 

- 

- 

17 

- 

- 

- 

(Zona  4) 

27.0 

- 

- 

4 

4 

- 

- 

90.0 

- 

• 

1 

1 _ : 

1 

- 

- 

8.3 

0.398 

0.9 

1 

- 

- 

- 

. 

4.3 

1  17.3 

- 

- 

3.18 

- 

- 

(Zone  5) 

27.0 

- 

18.3 

.  - 

- 

1.59 

- 

90.0 

1 

i 

20.0 

Hi 

• 

- 

• 

13.2 

0.542 

0.9 

! 

- 

0 

- 

- 

- 

4.5 

- 

14.3 

- 

0 

- 

- 

- 

(Zone  6) 

27,0 

- 

- 

■ 

3 

- 

- 

- 

90.0 

- 

18.7 

0 

“ 

* 

- 

Zone  4 

Zone  5 

Zone  6 

Huabsr 

of  Hazardous  Otbrls 

8.18 

1.59 

0 

Debris 

Dtiwlty  (4) 

'  98.2 

19.1 

0 

(1)  Ooly  the  low  launch  cnglca  w«r«  consldored  In  this  exaapla.  la  accual 
application  the  full  range  of  plaustbie  initial  launch  angles  should  be 
considered. 

(2)  In  this  illustrative  exaaiple,  the  debris,  kinetic  energy  and  height  were 
checked  only  at  cbe  aldspan  of  each  recovery  area.  In  practice,  the 
debris  energy  and  height  should  be  checked  et  all  positions  over  Che 
recovery  zone.  In  effect  then,  the  nethod  calls  for  checking  the  kinetic 
energy  for  every  particle  which  has  a  trajectory  which  intercepts  a  voluae 
3  a  (9.8  ft)  high  and  encompassing  the  entire  recovery  zone. 

(3)  velocity  is  given  if  the  debris  missile  did  not  pass  through  the  vertical  recovery  plane. 

*  2  2 

(4)  !fuaber  of  debris  aleoents  divided  bv  4.65  s  (30  ft  )  vertical  racovary  araa  and  normalized 
to  35.7  m^  (600  ft^). 
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CHAPTER  4 

OEKRIS  HAZARDS  METHODOLOGY  DEVELOPMENT  PLAN 


When  an  explosion  occurs  within  a  building,  structural  collapse  can  occur 
and  hazardous  debris  can  be  hurled  towards  other  occupied  areas.  Tl>e  mechanisms 
which  control  the  debris  formation  are  complex  and  vary  widely  for  different 
types  of  structures  and  different  parameters  associated  with  the  explosive 
contained  inside.  In  order  to  approach  the  problem  of  deriving  prediction 
methodologies  for  debris  hazards,  the  first  step  must  be  to  identify  the 
important  physical  parameters  which  have  an  effect  on  debris  characteristics. 
Identification  of  these  parameters  will  allow  one  to  organize  the  steps  to  be 
taken  in  an  investigation  of  debris  hazards. 

Debris  characteristics  are  defined  in  this  report  as  the  number  of  debris 
fragments,  debris  mass,  shape,  size,  velocity,  launch  angle,  range,  and  spray 
angle.  These  debris  characteristics  are  controlleJ  by  the  structural  details  of 
the  building  and  the  applied  internal  loads.  The  structural  details  depend  upon 
the  type  of  building  considered  and  from  what  portion  of  the  structure  the 
debris  originates.  The  internal  loads  are  dependent  upon  the  charge  size, 
shape,  location,  and  orientation  within  the  structure.  Also  the  effect  of 
charge  casing,  when  present,  can  affect  the  internal  loading. 

STRUCTURAL  PARAMETERS 

Debris  emanating  from  a  building  due  to  an  accidental  internal  explosion  is 
controlled  by  a  number  of  parameters  related  to  the  structural  details  of  the 
building.  Various  types  of  buildings*  are  used  to  house  ordnance  items 
including  reinforced  concrete  (laced  and  unlaced),  concrete  masonry 
construction,  and  steel  construction.  These  different  types  of  buildings  would 
each  present  different  classes  of  hazards  depending  on  the  failure  pattern.  This 
wide  range  of  construction  types  would  necessitate  a  very  broad  research  effort 
to  identify  debris  hazard  methodologies.  In  order  to  maintain  a  smaller 
research  effort  and  to  concentrate  on  the  type  of  construction  which  is  expected 
to  be  the  most  prevalent  for  ordnance  handling,  a  generic  debris  donor  structure 
should  be  selected.  This  report  will  concentrate  on  reinforced  concrete 
structures.  This  is  not  to  say  that  hazards  do  not  exist  for  structures  other 
than  reinforced  concrete. 


*The  contributions  from  equipment  and  furnishings  inside  the  building  are  not 
considered.  ' 
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Past  experience  has  shown  that  the  main  debris  hazard  from  explosions 
inside  reinforced  concrete  buildings  originates  frcm  the  walls.  This  may  not  be 
true  of  all  building  shapes;  however,  many  buildings  will  have  flat  roofs  from 
which  the  debris  is  projected  upward.  The  debris  will  fall  on  or  very  near  the 
original  site.  In  contrast,  wall  debris  will  be  projected  large  distances.  In 
addition,  the  wall  debris  emanates  primarily  in  a  airection  normal  to  the 
original  wall  surface.  This  is  shown  in  Figure  11.  By  concentrating  on  wall 
debris,  the  number  of  permutations  of  possible  debris  conditions  to  be  studied 
is  reduced. 

A  reinforced  concrete  wall  can  be  either  of  laced  or  unlaced  construction. 
Laced  construction  is  typically  used  in  the  design  of  a  structure  to  withstand 
the  applied  blast  loading.  The  specific  intent  in  the  design  of  a  laced 
reinforced  concrete  wall  is  to  control  the  collapse  and  breakup  of  the  wall. 

Some  debris  may  occur  in  the  form  of  backface  spall  or  scabbing.  For  occupied 
areas  exposed  to  these  debris  sources,  it  is  standard  practice  to  provide 
protection  in  the  form  of  spall  plates.  For  these  reasons  a  laced  reinforced 
wall  shall  not  be  considered  in  this  report.  Any  further  reference  to 
reinforced  concrete  will  pertain  to  unlaced  construction  unless  otherwise  noted. 

The  following  is  a  description  of  the  type  of  structure  considered  in  this 
report.  Typically,  the  wall  is  rigidly  attached  on  two  or  more  edges  with  equal 
reinforcing  running  in  both  the  horizontal  and  vertical  directions  and  on  both 
faces  of  the  wall.  The  concrete  cover  over  the  rebar  conforms  to  American 
Concrete  Institute  (ACI)  requirements,  0.019  meters  (0.75  in.)  for  inside 
surfaces  and  0.038  meters  (1.5  in.)  for  exterior  surfaces. 

Other  structural  conditions  could  have  an  effect  on  reinforced  concrete 
debris  formation.  Wall  reinforcement  can  be  tied  into  the  floor,  roof,  or  walls 
to  lorm  various  boundary  conditions  (fixed,  simple,  or  free).  The  type  of  wall 
boundary  conditions  will  control  the  yield  line  formation  which,  in  turn, 
affects  debris  formation.  The  amount  of  tie  in  rebar  and  wall  thickness  at  the 
perimeter  will  determine  whether  the  wall  fails  due  to  support  reactions  or  in 
bending.  Also,  stirrup  design  will  govern  whether  or  not  the  main  flexural 
rebar  is  allowed  to  develop  its  full  bending  resistance  without  premature 
failure  in  shear.  The  differences  in  types  of  failure  (support  reaction,  shear, 
or  bending)  could  cause  differences  in  debris  characteristics.  As  will  be 
suggested  in  Chapter  5,  the  model  test  series  will  include  a  test  fixture  with  a 
support  frame  in  which  the  debris  forming  wall  will  be  held.  The  number  of 
sides  of  the  wall  to  be  supported  in  the  frame  can  be  varied.  For  the 
investigation  suggested  in  this  report  only  fixed  boundary  conditions  should  be 
considered  at  this  time.  As  discussed,  shear  reinforcement  could  have  an  effect 
on  debris  formation.  Incorporating  shear  rebar  into  the  model  panels  to  be 
discussed  in  Chapter  5  would  be  costly.  Hence,  it  is  suggested  that  shear 
design  not  be  considered  for  the  initial  investigation  discussed  in  this 
report.  If  a  more  detailed  test  series  is  considered  at  a  later  date,  then 
shear  design  can  be  investigated  along  with  other  parameters  not  included  in  the 
initial  test  series.  It  should  be  noted  that  if  full-scale  tests  were 
conducted,  incorporation  of  shear  reinforcing  into  the  test  panels  would  also  be 
costly  and  a  scale  model  panel  including  shear  rebar  would  still  be  less 
expensive. 
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FIGURE  11.  DRAWING  OF  MEASURED  DEBRIS  DENSITY  OF  BUILDING  AND 
CRATER  MATERIAL  DURING  MODEL  TESTS‘^“ 
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For  the  reinforced  concrete  wall  described  above,  the  parameters  tf,  on 
considered  further  which  can  be  identified  as  controlling  the  debris 
characteristics  include  the  following: 


o 

o 

o 

o 

o 

o 

o 

o 

o 


Concrete  compressive  strength 
Aggregate  size 
Concrete  density 

Wall  dimensions  (height,  width,  thickness) 

Rebar  diameter 

Rebar  spacing 

Rebar /concrete  adhesion 

Rebar  yield  and  ultimate  strength 

Boundary  conditions,  i.e.,  wall  restraints 


These  parameters  and  their  relationship  between  each  other  and  debris 
characteristics  are  discussed  in  the  section  entitled  "Parameter  Relationship." 


BLAST  AND  FRAGMENT  LOAD  PARAMETERS 

The  formation  of  reinforced  concrete  debris  due  to  an  internal  explosion 
is  affected  by  the  applied  loading  which  is  a  function  of  charge  size, 
location,  and  the  presence  of  casing.  Shock  strength  increases  with 
decreasing  charge  standoff  and  increasing  charge  size.  Charges  close  to  the 
wall  but  off-center  in  the  room  will  have  a  shock' reverberat ion  pattern  which 
will  be  dominated  by  the  first  reflection.  Charges  centered  in  the  room  will 
have  secondary  reflections  which  are  also  important.  The  charge  size  and 
importance  of  secondary  shock  reflections  can  affect  whether  the  structure 
will  respond  impulsively,  dynamically,  quasi-statically,  or  under 
reverberation.  The  charge  size  and  location  within  the  structure  are  hence 
identified  as  important  parameters  controlling  debris  formation.  Room  volume 
affects  the  quasi-static  phase  of  an  explosion  and  is  also  considered  an 
important  blast  parameter.  Vent  area  is  not  considered  here  because  this  is 
usually  compared  to  wall  area,  and  wall  collapse  should  occur  before 
substantial  venting  for  large  charge  weights. 

The  shape  of  the  explosive  will  have  an  effect  on  the  forcing  function. 
High  explosives  can  be  in  boxes,  cylindrical  containers,  or  in  stacks  of 
munitions  of  a  variety  of  shapes.  Except  for  munitions,  nonspherical 
containers  are  expected  to  be  compact  (i.e.,  L/D  approximately  1)  for  ease  of 
handling.  Munitions  stacked  in  an  array  are  expected  to  be  in  a  compact 
arrangement.  To  consider  all  possible  permutations  of  charge  shape  would  be 
expensive,  particularly  when  a  spherical  or  compact  cylindrical  (L/D  =  1) 
shape  can  reasonably  represent  many  shapes  typically  encountered. 

Whether  a  charge  is  cased  or  uncased  will  have  an  effect  on  the  loading 
applied  to  Che  wall  surface.  The  casing  will  affect  the  blast  loading  and  can 
become  a  source  of  fragment  impact  on  the  wall  surface.  To  limit  the  type  of 
casing  to  be  studied  from  the  variety  of  casings  available,  it  is  suggested 
Chat  only  steel  casing  should  be  considered. 
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The  density  of  the  explosive  material  can  have  an  effect  i  .  the  loading  on 
the  wall  surface.  If  the  same  mass  of  explosive  is  packaged  iu  bulk  form,  a 
different  blast  field  is  expected,  given  an  explosion,  as  compared .with  that 
in  pressed  form.  Differences  are  also  expected  for  explosives  which  are 
packaged  as  individual  quantities  in  a  stack  (although  each  is  in  pressed 
form)  compared  with  a  single  pressed  unit  of  the  same  overall  shape.  The 
worst  case  loading  will  be  expected  from  the  compact,  pressed  explosive  and 
will  be  considered  under  this  report. 

In  summary,  the  explosive  parameters  to  be  considered  further  which  are 
identified  as  important  in  debris  formation  of  reinforced  concrete  walls  are: 

o  Explosive  mass 

o  Explosive  location  within  the  enclosed  structure 
o  Room  volume  (affects  quasi-static  phase) 
o  Explosive  shape  (only  compact  shape  will  be. considered) 
o  Casing  effects  (shock  interaction  and  jetting  effects  of 
warheads) 

o  Explosive  density  (only  pressed  explosive,  will  be  considered) 


PARAMETER  RELATIONSHIP 

The  parameters  previously  identified  as  being  important  to  debris 
formation  from  reinforced  concrete  walls  can  be  related  to  debris 
characteristics  (mass,  shape,  size,  velocity,  trajectory,  quantity,  range,  and 
distribution)  through  the  use  of  similarity  methods.  Similarity  is  a  powerful 
tool  for  organizing  parameters  into  terms  (called  pi  terms)  which  are 
physically  significant  to  the  problem  under  study.  Tlie  validity  and 
usefulness  of  similitude  analysis  is  well  documented  in  the  literature  and  has 
been  applied  to  a  wide  range  of  physical  problems  including  those  similar  to 
this  debris  problem  (Reference  34  includes  an  excellent  overview  of  similitude 
analysis).  Under  a  previous  study,  a  similitude  analysis  for  explosive 
fracturing  of  a  dividing  wall  was  developed. ^6  xhe  parameters  considered  in 
this  report  are  similar  to  those  studied  in  Reference  26,  however,  enough 
differences  exist  that  a  separate  similitude  analysis  was  necessary.  The 
results  of  Reference  26,  however,  were  used  as  a  guideline.  Similitude  (or 
model)  analysis  is  typically  used  in  conjunction  with  model  or  scaled  tests. 

A  model  test  series  will  be  suggested  and  discussed  in  Chapter  5  for  a  deb'.is 
hazard  study. 

The  pi  terms  derived  using  the  parameters  identified  in  the  development  of 
the  model  analysis  and  the  debris  characteristics  to  be  measured  are  given  in 
Tables  18,  19,  and  20.  The  tables  are  arranged  according  to  parameters  to  be 
varied  in  the  model  tests,  those  to  remain  constant,  and  the  parameters  to  be 
measured.  The  reason  for  varying  some  and  not  others  is  discussed  in  Chapter 
4.  Tables  18  through  20  include  important  parameters,  symbols,  and  dimensions 
(force  F,  time  T,  and  length  L),  Below  the  tables  are  the  nondimensional  pi 
terms  corresponding  to  the  above  parameters.  The  nondimensional  terms  were 
organized  to  reflect  physically  significant  relationships  thought  to  be  of 
importance.  Concrete  compressive  strength  is  comSined  with  wall  dimensions 
(which  govern  the  amount  of  concrete  in  compression  during  bending)  to  relate 
a  strain  energy  to  charge  energy.  The  yield  and  ultimate  strengths  of  the 


67 


NSWC  TR  85-114 


TABLE  18.  PARAMETERS  TO  BE  VARIED  DURING  MODEL  TESTING 


Parameter 


Symbol 


Dimension 


Wall  Thickness  X 

Rebar  Diameter  D 

Rebar  Spacing  S 

(horizontal  *  vertical) 

Charge  Energy  W 

Charge  Location  In  Room 

height  above  floor  h 

distance  from  wall  r 

distance  from  nearest  sidewall  l 

Charge  Casing  Thickness  d 


L 

L 

L 

FL 

L 

L 

L 

L 


a  XLH 


112  “  D/X 

113  «  s/x 

n4  -  h/H 

ils  =•  r/R 

Hs  -  l/L 

d^ac 


Note;  L,  H,  and  R  are  defined 
in  Table  19 
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TABLE  19.  PARAMETERS  TO 

BE  HELD  CONSTANT 

DURING  MODEL  TESTING 

Parameter 

Symbol 

» 

Dimension 

'  Room  Size 

height 

H 

L 

distance  £rom  debris  forming 

R 

L 

wall  to  back  wall 

distance  from  side  wall 

L 

L 

to  side  wall 

Aggregate  Size 

a 

L 

Concrete  Density 

FT^/L^ 

Concrete  Compressive  Strength 

a 

c 

F/l2 

Rebar  Density 

FT^/L^ 

Rebar  Yield  Strength 

•^y 

F/L^ 

Rebar  Ultimate  Strength 

0 

u 

F/L^ 

Gravity 

g 

L/T^ 

Standard  Air  Pressure 

Po 

F/L^ 

Strain  Rate 

• 

€ 

L/L/T 

na 

- 

(where  V  ■  HRL  •  room  volume) 

ng 

8K 

m 

K  MX 

His  -  ^ — 

Hu 

- 

a/x 

Hie  - 

w 

ni2 

> 

c  s 

ni7  =  et 

m 

(?  DS  (LH)^^^ 

y 

Note:  M  Is  defined  in  Tabli 

W 
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TABLE  20.  MEASURED  PARAMETERS  OR  DEBRIS  CHARACTERISTICS 


Symbol 


Launch  Angle  6 

Spray  Angle  ^ 

Debris  Velocity  v 

Debris  Mass  M 

Debris  Size 

length 

width  02 

thickness  03 

Debris  Range  Z 

Debris  Quantity  M 


Dimension 

L/T 

FT^/L 

L 

L 

L 

L 


ni8  *  9 


^21 

^22 

^^23 

^24 


H 

‘’c  XS 

(“1«2®3)  Pc 
M 


N 
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rebar  are  treated  similarly.  Rebar  spacing,  rebar  diameter,  and  aggregate 
size  are  combined  with  wall  thickness  which  for  these  parameters  is  considered 
to  be  the  most,  important  length  term.  Charge  location  within  the  room  is 
combined  with  the  corresponding  major  room  dimension  (i.e.,  charge  height 
above  floor  room  height).  Room  volume,  charge  weight,  and  atmospheric 
pressure  are  combined  relating  quasi-static  pressure.  The  parameters  in  Table 
20  are  also  arranged  with  similar  parameters  to  form  pi  terms.  Velocity  is 
squared  and  combined  with  mass  to  form  a  kinetic  energy  which  is  then  ratioed 
with  the  charge  energy.  Mass  is  combined  with  concrete  density  and  debris 
mass  (which  thereby  relates  debris  size  back  to  wall  dimensions).  Debris 
range  is  related  to  gravity  and  velocity  squared. 

All  of  the  pi  terms  discussed  above  are  suggested.  Rearrangement  of  terms 
can  be  made  if  more  suitable  relationships  are  found. 


SCALING  LAWS  FOR  REINFORCED  CONCRETE  MODEL 


Using  the  method  of  similarity,  it  is  possible  to  derive  the  scale  factors 
for  all  parameters  in  a  replica  model.  A  replica  Model  is  one  in  which  both 
the  prototype  and  the  model  structure  are  made  from  'he  same  materials.  In  a 
replica  model  all  geometric  dimensions  (lengths,  ti^  okness,  etc.)  are  scaled 
down  by  the  same  factor  X.  All  angles,  densities,  strengths,  strains,  and 
the  number  of  debris/fragmeiita  generated  will  be  invariant  or  the  same  as  in 
full  scale.  Time  and  impulse  in  the  model  will  scale  with  the  geometric  scale 
factor  X.  Mass  of  debris  and  energy  in  the  explosive  charge  will  scale  as 
X^.  Strain  rate  and  acceleration  due  to  gravity  will  scale  as  X“^. 

The  scale  factors  for  the  explosive  loading  and  subsequent  breakup  of  the 
reinforced  concrete  wall  are  summarized  in  Table  21.  All  parameters  listed  in 
Tables  18  to  20  will  be  scaled  properly  except  acceleration  due  to  gravity, 
debris  range,  and  strain  rate.  The  failure  to  scale  gravity,  which  should  be 
X  times  greater  in  the  model,  is  important  only  in  the  trajectory  of  the 
debris  generated  by  the  explosion.  Since  the  initial  launch  angle,  initial 
velocity,  and  debris  mass  are  properly  scaled,  it  should  be  possible  to 
calculate  the  equivalent  full-scale  ranges  using  these  initial  conditions 
measured  in  model  scale  experiments.  The  failure  to  scale  strain  rate 
properly  arises  because  strain  is  invariant  and  time  scales  as  X .  Thus 
strain  rates  in  the  model  will  be  X  times  larger  than  those  in  the 
prototype.  This  is  not  as  catastrophic  as  it  seems  because  steel  and  concrete 
are  not  strongly  sensitive  to  strain  rate.  The  ratio  of  "dynamic  strength"  to 
static  strength  for  these  materials  is  roughly  linear  with  respect  to  the 
logarithm  of  the  strain  rate.  In  an  explosion,  strain  rates  in  the  steel 
elements  and  in  the  concrete  may  vary  from  10”*^  to  lOO/sec.  Over  this  range 
those  materials  may  have  an  increase  in  "dynamic  strength"  from  1.1  to  1.8 
times  the  static  strength.  In  a  1/lCth  scale  model  the  strain  rates  in  a 
similar  explosion  may  be  10®  to  1000/sec.  The  "dynamic  strength"  increase 
in  the  model  may  then  be  on  the  order  of  1.2  to  2.4  times  the  static 
strength.  From  this  discussion  it  can  be  seen  that  strain  rate  effects 
between  the  model  and  prototype  structure  are  real;  however,  these  effects 
probably  will  be  difficult  to  discern  given  the  usual  scatter  inherent  in  this 
type  of  tfe  •^ing. 
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TABLE  21.  MODEL  LAW  FOR  DIVIDING  WALL  FRAGMENTATION 


Parameter 

Replica  Scaling  Law 

Lengths 

X 

Angles 

1.0 

Densities 

1.0 

Strengths,  moduli 

1.0 

Poisson's  ratio 

1.0 

Strains 

1.0 

Strain  Rate 

x-1 

Velocities 

1.0 

Acceleration  (gravity) 

x-^ 

Mass 

x^ 

Reinforcement  ratio 

1.0 

Explosive  energy 

x3 

Pressure 

1.0 

Impulse 

X 

Time 

X 

Number  of  fragments 

1.0 
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NUMERICAL  RANGE  OF  PARAMETERS  TO  BE  INVESTIGATED 

As  discussed  in  the  previous  section,  it  is  suggested  that  the  parameters 
listed  in  Table  18  be  varied  during  model  testing.  The  parameters  in  Table  19 
will  be  held  constant  £or  all  tests.  The  reason  these  parameters  are  held 
constant  is  because  the  effect  of  their  variation  on  debris  formation  for 
structures  of  interest  is  expected  to  be  less  significant  (except  for  room 
dimensions)  than  parameters  listed  in  Table  18.  The  parameters  in  Table  19 
are  also  somewhat  constant  in  typical  construction,  i.e.,  aggregate  size, 
concrete  density,  and  rebar  strengths  do  not  vary  greatly  in  common 
construction.  The  number  of  permutations  can  be  greatly  reduced  by  holding 
these  parameters  constant.  Consideration  of  a  constant  room  size  is  strictly 
an  economical  decision.  The  model  test  setup  to  be  discussed  in  Chapter  5  is 
envisioned  to  be  a  reusable  structure  (except  for  the  debris  forming  wall), 
and  it  is  economically  desirable  to  maintain  a  constant  room  size. 

The  range  over  which  the  parameters  in  Table  18  are  to  be  varied  was 
chosen  to  reflect  that  expected  to  be  found  in  existing  structures.  The 
ranges  are  given  in  Table  22.  A  survey  of  existing  structures  was  not  made 
for  use  in  determining  the  range  of  variance  in  the  table,  but,  these  values 
were  chosen  for  what  past  experience  has  shown  to  be  typical.  Any 
disagreement  with  these  values  can  be  incorporated  by  expanding  or  contracting, 
the  range  of  parameters  to  be  considered  as  listed  in  the  table 
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TABLE  22.  FULL-SCALE  RANGE  OF  VARIATION  TO  BE  APPLIED  TO  PARAI<1ETERS 
IN  TABLE  18  DURING  MODEL  TESTING 


Wall  thickness 

Explpslve  mass 

Casing  (one  thickness) 

Charge  Position  in  room 

height  above  floor 
distance  from  wall 
distance  from  sidewall 

Rebar  size 

Rebar  spacing 


0.3  -  0.46  meters  (12  -  18  in.) 
150,-  4500  kg  (330  -  10,000  Ib^^) 
Bare  or  cased 

0  -  1/3  room  height 
varied 
centered 

//5  -  //8 

0.15  -  0.46  m  (6  -  18  in.) 
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CHAPTER  5 

recommf.nded  model  tests 


To  better  define  the  debris  characteristics  of  reinforced  concret  ; 
buildings  and  walls  exposed  to  blast  loading,  it  is  recommended  that  a  small 
scale  study  be  initiated.  The  use  of  model  tests  can  provide  a  large  data  base 
at  a  reduced  cost  compared  to  full-scale  testing.  With  the  use  of  the 
similitude  analysis  discussed  in  Chapter  4  the  effect  of  the  various  parameters 
on  debris  formation  can  be  studied  using  model  test  data.  Model  testing  is 
widely  used  for  problems  similar  to  debris  hazards.  Modeling  is  well  understood 
and  documented  in  the  literature  (References  35  and  36)  where  the  validity  of 
properly  managed  model  tests  are  unquestioned.  The  model  used  for  testing 
should  simulate  the  prototype  as  closely  as  possible  for  all  important  physical 
parameters.  These  conditions  are  easily  satisfied  except  for  the  gravity  term 
which  unfortunately  does  not  scale.  This  handicap  only  affects  the  range  of  the 
debris,  i.e.',  the  range  term  does  not  scale  properly.  All  other  measured 
parameters,  however,  do  scale.  During  model  testing  several  parameters  vill  be 
varied.  The  overall  test  program  can  be  performed  as  k  series  of  small  cest 
programs  where  one  or  two  ^f  the  parameters  will  be  varied  with  the  remaining 
parameters  held  constant.  At  time  of  completion  of  the  first  series  of  tests, 
another  series  would  be  initiated  wiph  a  different  set  of  parameters  varied. 

For  validation  of  the  small  scale  tests,  several  model  scales  including  full 
scale  should  be  investigated.  Table  23  outlines  a  series  of  tests  designed  to 
provide  accurate  and  reliable  fragment/debris  data  for  reinforced  concrete 
structures  subjected  to  internal  loads.  As  mentioned  earlier,  when  the  final 
s'-’ection  of  the  donor  system  is  made,  the  test  plan  and  corresponding  modal 
test  series  can  be  organized  in  a  manner  similar  to  that  discussed  here. 


SCALE  MODEL  ENCLOSURE 

It  is  suggested  that  the  tests  to  be  performed  on  reinforced  concrete  w^lls 
make  use  of  a  scale-model  enclosure.  This  enclosure  would  consist  of  a  reusable 
corxrete  base  slab  (to  which  a  frame  or  fixture  for  use  in  supporting  the  roof 
and  wall  3  would  be  attached),  a  blowout  roof,  and  three  frangible  walls'.  The 
fourth  wall  will  be  the  test  wall.  The  test  wall  will  be  made  of’  reinforced 
concrete  and  would  be  supported  on  either  one  or  three  sides  by  the  reusable 
support  frame. 

SwRI  fabricated  1/6  scale  model  reinforced  concrete  walls  for  a  program 
funded  by  aRRADCOM,  entitled  "Explosive  Fragmentation  of  Dividing  Walls. 

In  fabricating  the  model  walls,  14  gauge  steel  wire  was  used  to  simulate  the 
reinforcing  bars.  The  wire  was  chosen  to  have  equivalent  strength  properties  of 
the  rebar  being  modeled,  and  the  gauge  of  the  wire  was  a  scaled  equivalent  of 
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the  rebar  size.  Because  rebar  is  not  smooth  and  bonding  to  concrete  is 
importan'.,  the  wire  was  crimped  to  incorporate  this  feature  ip  the,  model.  Rebar 
was  located' in  both  the  horizontal  and  vertical  direction  without  lacing.  Molds 
Cor  the  concrete  walls  were  fabricated  into  a  rectangular  plywood  frame  which 
was  designed  to  be  reusable.  The  frame  had  a  series  of  holes  on  the  side  which 
held  the  rebar  in  its  proper  place.  Concrete  was  poured  into  the  model  such 
that  a  lav»r  of  concrete  covered  the  rebar  on  both  faces.  SwRI  recommends  this 
type  of  fabrication  process  for  use  in  the  proposed  program  if  small  scale 
models  are  used.  For  larger  scale  models  .actual  rebar  is  available  for  use  in 
fabrication  of  the  test  wall.  The  smallest  rebar  available  commercially  is  No. 

2  bai ;  however,  it  is  only  obtainable  in  plain  round.  Tl\is  limits  the  scale  of 
the  model  if  actual  rebar  is  to  be  used.  * 

The  frangible  walls  and  roof  can  be  fabricated  of  light  yet  sturdy  material 
such  as  plywood  or  pressboard.  The  frangible  panels  would  be  supported  in  such 
a  way  that  they  would  remain  in  place  long  enoug’u  to  reflect  the  shock  waves 
generated  by  the  explosion.  This  would  subject  the  test  wall  to  a  realistic 
forcing  function.  The  frangible  panels  can  be  covered  by  dirt  to  simulate  a 
full-scale  mass  if  uncovered  frangible  wall  response  is  determined  to  be  too 
rapid.  The  dir.t  can  be  heaped  against  and  on  the  model  or  poured  between  two 
panels  of  proper  spacing  to  obtain  the  correct  scaled  mass.  When  the  frangible 
walls  and  roof  (either  dirt  covered  or  uncovered)  do  fail,  the  fact  that  they 
are  made  of  light  materials  will  eliminate  hazardous  debris. 

The  test  wall  should  be  color  coded  on  both  the  front  and  back  face.  It  is 
suggested  that  at  least  four  colors  be  used  on  each  face.  To  identify  sectidns 
of  the  test  wall,  small  number  labels  should  be  glued  on  the  wall  in  an  array 
within  each  colored  area.  Previous  model  teats  have  incorporated  this  color 
coding  and  labeling  scheme  with  success. The  tests  should  be  made  with 
the  wall  facing  a  debris  recovery  pit  which  is  a  sand  bed  large  enough  to 
receive  all  expected  debris. 


INSTRUMENTATION 

It  is  suggested  that  both  real  time  and  high-speed  cameras  be  used  to 
record  the  event.  The  cameras  can  be  placed  at  tlie  model  end  of  the  recovery 
pit  in  order  to  record  the  debris  immediately  after  the  explosion.  A  ruled 
background  panel  can  be  used  to  serve  as  a  backdrop  for  the  car;era  recordings. 
Ruled  lines  on  the  panel  will  give  a  reference  frame  iii  the  film.  This  type  of 
setup  was  used  with  the  model  tests  described  in  Reference  26.  Figure  12  is  an 
example  of  this  kind  of  test  setup. 

Pressure  gauges  can  be  located  on  the  test  enclosure  floor  if  a  record  of 
the  blast  environment  is  desired.  Because  the  test  wall  will  collapse,  no 
pressure  history  recordings  can  be  made  of  tlie  wall  surface  itself.  Blast 
measurements  can  be  made  if  a  steel  panel,  designed  not  to  fail,  is  placed  in 
the  test  frame  with  an  attached  array  of  blast  gauges.  This  type  of  test  setup 
is  referred  tq  as  a  loads  model.  A  minimum  of  five  gauges  evenly  spaced  on  the 
plate  is  suggested.  In  addition  to  those  on  the  plate,  blast  gauges  should  be 
located  on  the  floor  (in  the  same  location  as  for  debris  model)  for  use  as 
reference  gauges.  The  loads  model  would  require  a  series  of  tests  which  include 
all  charge  sizes  and  locations.  A  comparison  of  the  pressure  histories  from 
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reference  gauges  (on  Che  floor)  could  be  nade  between  a  corresponding  loads  test 
and  debris  test.  This  would  allow  deCermiTiaC'ion  of  whether  the  applied  loads  to 
Che  wall  are  similar. 


DEBRIS  DATA  TO  BE  RECORDED 

Measurements  made  for  each  debris  cest  i^hould  include  the  recorded  values 
corresponding  to  the  parametets  identified  in  Table  21  for  a  reinforced  concrete 
wall.  Other  debris  donors  would  have  similar  characteristics  to  be  measured. 
Debris  mass  and  dimensions  should  be  measured.  The  presented  area  can  be 
estimated  in  two  ways.  The  best  method  consists  of  an  icosahedron  gauge  which 
uses  phoCodetectors  to  measure  the  "shadow”  of  the  debris  in  each  of  16 
different  orientations.  The  16  different  readings  are  averaged  to  yield  an 
"effective  presented  area."  Unfortunately  this  device  is  quite  expensive  and 
only  two  such  units  are  in  existence.  A  simple  and  cheaper  method  consists  of 
measuring  three  orthogonal  dimensions  (L,  W,  T)  of  Che  debris.  From  these 
readings  an  "effective  presented  area"  is  estiraated  as  (L  x  W  x  T)^/^. 

Debris  range  should  be  measured  and  recorded.  An  attempt  should  be  made  to 
locate  Che  origin  of  Che  debris  on  the  wall  through  use  of  the  color  coding  and 
numbering  system.  Debris  trajectory  and  veloci*'y  may  be  obtained  from 
high-speed  film  recordings  and  perhaps  flash  techniques  or  breakwire  systems. 

The  number  of  debris  particles  per  test  should  be  recorded.  If  very  small 
debris  makes  this  difficult,  Chen  debris  above  a  certain  size  can  be  counted. 

The  number  of  debris  elements  recovered  versus  distance  from  Che  test  wall  can 
be  measured.  The  number  of  debris  particles  versus  debris  mass  can  be  counted. 
These  parameters  can  be  obtained  from  what  is  corimonly  called  a  missile  map 
which  includes  debris  location,  total  number  of  debris,  and  debris  size. 


ADDITIONAL  TESTS  IN  SUPPORT  OF  MODEL  TESTS 

In  Chapter  4  Che  important  parameters  affecting  debris  characteristics  were 
discussed.  One  of  the  parameters  which  affect  debris  range  is  gravity  (ir23. 
Table  20),  which  is  not  properly  scaled.  It.  order  tn  determine  a  full  scale 
range,  additional  analysis  of  the  test  data  is  required.  Using  Che  debris, 
velocity,  mass,  shape,  launch  angle,  and  spray  angle,  a  corresponding  range  can 
be  determined  using  existing  trajectory  .solutions.  The  only  missing  parameter 
is  a  proper  drag  coefficient  and  debris  presented  area. 

Since  there  is  no  method  for  predicting  the  drag  coefficient  or  the  drag 
area  for  tumbling  debris,  it  is  common  practice  to  assume  a  stable  orientation 
for  the  debris  missile,  calculate  the  drag  areas  tor  this  orientation,  and 
obtain  Cd  from  handbooks  or  standard  references.  In  addition,  values  for  Cj) 
are  given  for  only  a  limited  number  of  orientations,  i.e.,  edge-on,  side-on, 
etc. 37  Coefficients  for  other  debris  missile  orientations  such  as  10“ 
off-normal,  20“  off-normal,  etc.,  are  not  available.  In  reality,  debris 
missiles  will  have  any  of  numerous  orientations  including  tumbling.  Test  data 
covering  a  range  of  debris  orientations  would  be  very  useful  to  obtain  a 
prediction  model.  This  will  require  an  experimental  program  varying  the 
orientation  and  rate  of  tumble  of  typical  debris  fragments.  This  type  of 
program  could  be  performed  using  debris  missiles  with  known  dimensions  fired  out 
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of  an  air  gun  with  a  certain  orientation,  and  at  a  known  launch  angle  and 
initial  velocity.  The  range  traveled  by  the  missile  could  be  recorded  and 
compared  to  that  predicted  by  the  standard  trajectory  codes.  By  varying  the 
initial  debris  orientation,  ' the  resultant  debris  range  could  be  corrected  to 
account  for  the  drag  of  the  missile.  A  similar  analysis  could  be  performed  for 
tumbling  debris  by  introducing  a  tumbling  medium  and  then  changing  the  tumbling 
rate. 
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CHAPTER  6 

ANALYSIS  OF  TEST  DATA 


test  data  collected  as  discussed  in  Chapter  5  should  be  part  of  a 
parametric  study  aimed  at  developing  methodologies  for  the  reasonable  prediction 
of  hazards  associated  with  debris  from  blast  loaded  reinforced  concrete  walls. 
Proper  analytical  models  would  be  useful  as  input  to  the  NOHARM  program  for  use 
in  debris  hazard  predictions.  Generalization  of  experimental  results  is 
difficult  due  to  the  response  dependence  on  fabrication  and  explosive 
parameters.  To  date  there  are  no  established  prediction  equations  for 
calculating  the  effect  that  various  parameters  have  on  concrete  wall  break-up. 

In  Reference  26  a  model  analysis  was  used  to  design  the  scale  model  tests 
and  to  interpret  the  test  results.  Based  on  the  limited  tests  conducted  on  this 
program,  preliminary  prediction  models  were  formed  for  greatest  debris  velocity, 
longest  debris  range,  number  of  debris  elements,  and  largest  recovered  mass  as  a 
function  of  the  scaled  total  impulse  applied  to  the  wall.  Examples  of  the  types 
of  results  for  dividing  walls  supported  on  three  sides  are  given  in  Figures  13 
through  16.  Although  these  results  are  not  directly  applicable  to  the  situation 
of  a  naval  munition  exploding  within  a  structure,  these  results  are  indicative 
of  the  type  of  empirical  analysis  which  is  possible. 

Thus  far,  the  analysis  which  can  be  performed  on  maximum  responses  obtained 
in  the  tests  has  been  presented.  A  similar  analysis  can  be  performed  oi 
statistical  distributions  of  measured  data.  For  example,  consider  the  case  of 
the  mass  of  the  debris  fragments  collected  in  the  dividing  wall  tests. In 
these  experiments  all  significant  debris  generated  was  collected  and  weighed.  ' 
The  mass  distributions  were  prepared  by  choosing  mass  intervals  which  were 
approximately'  log  linearly  distributed.  The  number  of  debris  missiles  with  a 
mass  greater  than  the  specified  mass  intervals  was  counted.  Using  these  data, 
plots  similar  to  the  one  given  in  Figure  17  were  prepared.  The  data  from 
several  tests  were  plotted  in  Figure  18.  In  this  case,  the  various  curves 
represent  different  explosive  loadings  applied  to  the  wall.  Of  course,  the  top 
curve  represents  the  most  highly  loaded  wall,  and  the  bottom  curve  represents 
the  wall  with  the  lowest  loadings.  Each  of  the  curves  has  a  similar  shape  which 
suggests  that  it  should  be  possible,  through  proper  selection  of  dimensionless 
ratioSj,  to  correlate  empirically  all  of  the  mass  distributions  into  a  single 
curve.-*® 

The  above  are  suggestions  of  the  type  of  analysis  that  can  be  performed 
with  the  test  data.  Similar  analyses  are  expected  in  conjunction  with  the 
actual  test  series.  The  use  of  similarity  analysis  provides  a  powerful  tool  in 
development  of  predict! -n  methodologies.  The  data  accumulated  will  have  to  be 
reduced  for  use  in  NOHARM  and  should  be  formatted  as  required  by  the 
methodologies  described  in  Chapter  3. 
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FIGURE  15.  NUMBER  OF  DEBRIS  MISSILES  PRODUCED  AS  A  FUNCTION  OF  'UE 
IMPULSE  FACTOR  FOR  THREE-SIDE-SUPPORTED  PANET,s26 
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FIGURE  16.  LARGEST  RECOVERED  MASS  AS  A  FUNCTION  OF  THE  IMPULSE 
FACTOR  FOR  THREE-SIDE-SUPPORTED  PANELS26 
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